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1 Introduction. 



Although the experiments of this paper, statistical and numerical, were undertaken in pursuit 
of a goal not widely shared, they may be of general interest since they reveal a number of 
curious properties of the two-dimensional Ising model that had not been previously observed. 

The goal is not difficult to state. Although planar lattice models of statistical mechanics 
are in many respects well understood physically, their mathematical investigation lags far 
behind. Since these models are purely mathematical, this is regrettable. It seems to us that 
the problem is not simply to introduce mathematical standards into arguments otherwise well 
understood; rather the statistical-mechanical consequences of the notion of renormalization 
remain obscure. 

Our experiments were undertaken to support the view that the fixed point (or points) of 
the renormalization procedure can be realized as concrete mathematical objects and that a 
first step in any attempt to come to terms with renormalization is to understand what they 
are. We have resorted to numerical studies because a frontal mathematical attack without 
any clear notion of the possible conclusions has little chance of success. We are dealing with 
a domain in which the techniques remain to be developed. 

A fixed point is a point in a space of presumably an infinite number of dimensions; so 
this point and all other points of the space are defined by an infinite number of coordinates. 
Some will presumably be superfiuous, so that the total space is realized as a submanifold 
of some larger coordinate space. The total space will be the carrier, in some sense, of 
the renormalization transformation, but the transformation will not appear explicitly in 
this paper. The point does! The implicit condition on each quantity serving as a possible 
coordinate of the fixed point is that, at the fixed point itself, it remains invariant under 
renormalization and that, at a critical point of any model within the class considered, its 
value approaches a limit under repeated renormalization because renormalization drives the 
critical point to the fixed point. Since repeated renormalization is in coarsest terms nothing 
more than passage to larger and larger blocks or to smaller and smaller mesh, the condition 
is that the quantity has a meaning as the mesh length goes to zero, the dimensions otherwise 
remaining the same. For percolation this is a property of crossing probabilities. Our point 
of view is that any such quantity is a candidate as a possible coordinate in the space of 
the fixed point. Rather than a single numerical quantity we can consider several at once, 
which amounts in the customary mathematical way to considering objects lying in some 
given space, finite-dimensional or infinite-dimensional, for example, a space of probability 
distributions, and if these objects satisfy this criterion, thus if, for each model at the critical 
point, they tend to a limit as the mesh goes to zero, then this limit or rather its coordinates 
in the given space can also serve as coordinates of the fixed point. Such objects are described 
in the paper. 

There are at least two possibilities: one modeled on the considerations for the free boson 
of 0; the other on the crossing probabilities for percolation ||LPPS| , [LPS|| . The second 
possibility was suggested to us by Haru Pinson to whom we are grateful. Thus to each form 
M of the Ising model (taken at the critical temperature) we will attach two points pi:,{M) 
and p(7(M), each defined by an infinite number of coordinates. Both are, in so far as this 
can be confirmed by experiments, universal and conformally invariant in the sense of | LPS |. 
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It is unlikely that these two points are independent. One set of coordinates may well be 
deducible from another, but we have not examined this possibility. 

Crossing probabilities may or may not be peculiar to a few models. The evidence for 
their conformal invariance and universality is not difficult to present and appears in Section 
^. There is, however, one point to underline. The Ising model is considered in regions that 
may be bounded or unbounded. Crossing probabilities are defined for crossings within a 
region that may or may not coincide with the region in which the model is considered. It 
may be smaller. In contrast to crossing probabilities for percolation, in which there is no 
interaction, those for the Ising model depend on both the region in which the crossings are 
allowed to occur and the region of thermalization on which the Ising model is considered. 
Conformal invariance refers to the simultaneous action of a conformal mapping on the pair 
of regions. 

The coordinates modeled on the free boson should, on the other hand, be available for a 
large class of models. Their definition is, in principle, quite general, but we have confined 
ourselves to the Ising model. The states a of the Ising model take values in the set {±1} 
which is contained in the set of all complex numbers z with |z| = 1. This set in turn is covered 
by the line z = exp[27rix) . We simply develop the circle on the line. We first assign to one 
site Pq in the lattice the value /i(po) = and then choose for all other p the value h{p) = mn, 
m G Z, so that exp{ih{p) — ih{po)) = a{p)/a{po). Of course, there has to be more method 
than that. For example, for the square lattice we introduce clusters: maximal collections 
of lattice sites of the same spin that are connected through bonds joining nearest-neighbor 
sites. Each cluster is surrounded by a curve constructed from edges of the dual lattice and 
this curve separates it from all other clusters. To each of these curves an orientation is 
assigned randomly and, for nearest neighbors p and q, we set h{p) — h{q) = dzTr, the sign 
being determined by the relative orientation of the bond from ptoq and the curve it crosses. 
If it crosses no curve, then a{p) = a{q) and we take h{p) = h{q). Thus to every state a are 
attached several functions h, but h determines a up to sign. For a finite lattice the measure 
on Sj, the set of all possible h is taken to be such that the measure, mj, on assigns the 
same mass to all points lying above a given a. Their sum is one-half the mass of a. 

Fix now a bounded planar region D and consider the Ising model in this region with 
respect to a square lattice whose mesh a approaches 0. Since the model is to be critical, the 
contribution to the Boltzmann weight of a pair of neighboring spins is 

e-^'''''\ sinh(2J) = 1, J = .440687. 

Each h is in effect a function on the whole region if we take the value in the open square of 
side a about the site p to be h{p). (The ambiguity at the boundary is disregarded here; it has 
to be confronted in various ways from experiment to experiment.) If C is a (smooth) curve 
in R, which can lie entirely in the interior or run entirely or partially along the boundary, 
then we can restrict each function in Sj to C. This yields a set of points, each carrying a 
mass, in the set Vc of Schwartz distributions on C, and thus a probabihty measure on Vc- 
Experiments to be described in Sections ^ and ^ suggest strongly that this measure has a 
limit as the mesh tends to and that the limit is universal and conformally invariant. This 
is perhaps the most important conclusion of the paper. 



3 



These measures have surprised the authors more than once. When C is the boundary, 
the measure has a number of properties, to which we devote considerable attention, that 
suggest it is gaussian. It is not. 

There is no reason to restrict ourselves to planar regions and we begin our study with 
the cylinder, because the ambiguities at the boundary are then absent. A long cylinder 
(effectively semi-infinite) is, provided we stay close to the end, to be regarded as equivalent 
to a disk. The simplest conformally invariant distribution on the set of distributions on the 
boundary of a disk is the gaussian distribution with respect to the quadratic form defined 
by the Dirichlet form. For a function Lp this form is obtained by extending (y9 to a harmonic 
function (p in the interior and then taking 

The first experiments described in Section ^ strongly suggest that the measure given by our 
construction is in some respects very similar to this gaussian with a constant g = Qb that 
appears to be universal; the final experiments of that section show, however, that it differs 
in important respects from a gaussian. If C is interior, the measure on Vc is in no respect 
similar to a gaussian. 

Our construction is different from but not unrelated to familiar constructions relating the 
Ising model to SOS models. For the Ising model on a triangular lattice our construction is 
equivalent in many respects to the usual one for the 0(l)-model. In particular, it is expected 
that in the plane 

mx)-hmf) ~ -inb-oi. 

9 

The constant g = gj is expected to be 4/3, at least for the triangular lattice, but the two 
constants gj and gs are not equal. 

gi = 4/3, gs = 1.4710. 

Although gi is usually defined only for the triangular lattice, it can be defined in general. 
We suppose that it is universal, but we have not examined this carefully. The pertinent 
experiments are briefly discussed in Section |^. The conclusion, which will be reinforced 
more than once as we proceed, is that the comparison with the free boson undertaken in this 
paper is quite different than the usual one. 

There is no reason that the two constants gi and gs should be equal. One refers to interior 
behavior in the bulk, the other to behavior on the boundary. Moreover, as it turns out, they 
refer to different aspects of a construction that leads to nongaussian measures with some 
gaussian behavior. Although a departure from the conventional view, it could be argued (we 
do not attempt to do so here) that for questions of renormalization the constant gs may be 
every bit as important as gj, or, much better, that the distributions on curves of Section]^ are 
at least as important for renormalization as asymptotic behavior because renormalization, 
at least as it is often presented, entails the fusion of bounded regions along their boundaries 
(which may or may not partially coincide with that of the region of thermalization). The 

4 



measures on Vc were originally examined only for curves on or close to the boundary. They 
appear, somewhat to our surprise, to be of interest even in the absence of a boundary. Indeed 
it may turn out, with hindsight, that the numerical arguments towards the end of Paragraph 
3.2 for the existence of nontrivial and conformally invariant measures on interior curves are 
at least as important as the other results, argued more elaborately and with more detail, of 
Sections and 

Our point of view would not be at all persuasive if there were no sign in our fixed point p^:;, 
thus in the measures on Vc, of the critical indices 0, 1/2 and 1/16. It is seen in Paragraphs 
and B.2| that these measures do contain information about critical indices. Section HI, in 



which we describe another manifestation of the index 1/16 as well as an interpolation of a 
formula of Cardy, is also an essential part of the paper. 

The final section is less important. It contains a few observations that provide some 
perspective on the definitions of the paper. First of all, the construction of h is by no means 
canonical. There are alternative constructions decribed in Section |^. We can allow jumps 
other than ±7r, in particular several jumps nvr, n odd, with equal or different probabilities. 
They lead to different values of gi and to a measure on Vc with little resemblance to a 
gaussian. 

The possibility of not using clusters in our sense but the clusters of Fortuin and Kasteleyn 
that appear in the high-temperature expansion of the Ising (or more generally the Potts) 
model also suggests itself. Such clusters can also be used to define the crossing probabilities. 
They lead to different measures and to different crossing probabilities, whose universality 
and conformal invariance we have not tested. 

Finally we point out that the results for the distributions appear to remain valid at 
infinite temperature with, of course, a different value for the parameter appearing in the 
distribution. 

We are grateful to Michael Aizenman, Christian Mercat, Haru Pinson and Thomas 
Spencer for observations and suggestions that have been useful to us during the course of 
these investigations.^ Our most important debt is to the statistician Christian Leger whose 
counsel and advice were of enormous help in the preparation of Section |^. 

Questions and comments of the referee, to which in particular the discussion relating 
statistical dependence and the two-point correlation function in Paragraph is a response, 
led to at least one important modification of the conclusions in the original version of the 
paper. 



^RPL is grateful to Alois Scharf for making the bank of computers of the Mathematisches Institut der 
Universitat Bonn available to him at an important moment in the experimentation. 
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2 Distribution of h at the boundary. 



2.1 The free boson on domains with boundary conditions. 

The partition functions of a free boson 0, with compactification radius R, are famihar objects 
when the domain is a torus, or a rectangle with the field satisfying Dirichlet boundary 
conditions, thus equal to on the boundary. For a general Riemann surface with boundary 
and for an arbitrary specification of the field at the boundary, it may still be possible to 
describe the partition functions explicitly (see [ pG|| ). For a cylinder we use the formula of 



Q. As it suggested some of the statistical quantities for the Ising model studied in this 
paper, we review this formula. 

The cylinder is described as the quotient of the region < < — In g, 1 > g > 0, in the 
complex plane by the transformations generated hYZ-^z + uj,u! = 2i7r. If the fundamental 
domain is chosen to be < < — Ing and < < 2m, the map w identifies 
the cylinder with the annulus of outer radius 1 and inner radius q. We shall use freely the 
terminologies associated with the cylinder and with the annulus. Observe that q is close to 
zero for long cylinders. The angle 6 is used as the parameter on both the inner and the outer 
boundary. 

The extremal fields d) on the domain are real harmonic functions 



bo + a\n z + bin z + 



The boundary conditions fix the restriction of to the boundary. On the inner circle 
where z = ge*^ and z = ge~*^, this restriction is 

(j).^{e) = (j)o + {a + b)\nq + ie{a - b) + Y1 ^k^^"'"^ 

with the reality condition a^^ = af and on the outer circle 

^out{9) = ^o + ^0{a-b)+J2b^e''' 

with b^). = 6f . (The superscript stands for boson.) The compactification condition does not 
require to be periodic but imposes a milder condition: (j){e'^'^'^ z , e~^*'^z) = (j){z, z)—2'irnR, n G 
Z, thus {a — b) = inR, n E 1i. Since the Lagrangian function (||) does not depend on the 
term (po, this constant can be set to zero. Therefore only the difference of the constant terms 
in and (pont is of significance and we choose to parametrize it with a variable x G [0, 27tR) 
and an integer m G Z: 

-{a + b)\nq = X + 2TrmR. (2) 
The partition function on the cylinder with the boundary values of (p specified by 



6out, or equivalently by x, {af } and {&f }, is a product of three terms [L] 

^(0in,0out) = ^(a:,{af},{fef}) (3) 
= A-^Z,{x)Z,{{a^},{b^}) (4) 
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where A is the (^-regularization of the determinant of the Laplacian for the annulus. It is 
given by A = —iTrfij) where q = e*'^'^ and rj^r) = e™'^!^'^ 11^=1 (-'- ~" e^*'"'^'^) is the Dedekind 
ri function. Since this factor is independent of the boundary data, it will be disregarded. 
The crucial terms here are the two other factors ZAx) 



E 



q^^" ^ (5) 



(6) 



andZ2({af},{6f}) 

If measurements are made disregarding the variable x, only 

1-277 R 

Zix,{a^}Abk})dx 

is of importance and this gives, after proper normalization, a probabilistic measure on the 
space of boundary data parametrized by ({af }, })• The mixing of the boundary data at 
both extremities becomes more and more intricate when q approaches 1 or, in other words, 
when the cylinder becomes a narrow ring. When q is taken to zero, the measure simplifies 
as it becomes the product of two terms, each one depending on {af } or {6f }. Moreover, in 
the limit q = 0, the probabilistic interpretation of Z is simply that of the gaussian measure 
in the variables and b^. 

Even though the Coulomb gas provides a description of the minimal models, we do not 
know of any similar explicit formula for the partition functions of these models for general 
boundary conditions, although Cardy's paper ||C1|| treats explicitly the case of conformally 



invariant boundary conditions. There are indeed only a finite number of these, and one of 
the difficulties addressed in this paper is how to introduce continuously varying conditions. 
Nonetheless we proceed boldly using the partition function (^) as a guide for the Ising model. 
In contrast to the free-boson model, the Ising model defined on a graph Q does not have a 
field taking its values in the whole real line that we could easily identify with - the spin 
field a takes its values in {+1, — !}• Starting from the spin field a, defined on the sites of a 
(finite) graph Q, one can construct the function h as described in the introduction. It is such 
that, if p and q are joined, then h{p) — h{q) = ±n if a{p) ^ a{q) and h{p) = h{q) otherwise. 
If the graph Q is embedded in a surface D, for example, a cylinder or R^, this function h can 
be extended to a function locally constant on D except on the edges of the dual graph where 
it has jumps. The Ising measure m/ on the space of configurations on the graph Q of mesh 
a endows the (finite) set of possible functions h with a (discrete) probability measure. 
(As observed above, this measure is such that m'}^{h) = 2mj{a)/N^ where N^- is the number 
of distinct /I's that lead to a and — cr.) 

Take the graph to be the subset of the lattice aZ^ of mesh a = 1/LV formed by the 
points {am, an), < m < LH , Q < n < LV . We identify upper and lower edges and 
regard the graph as a subset of the cylinder: z = m + in ^ exY){—2Tiz/LV). How can we 
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compare to the field-theoretic measure of the free boson? Using the same letters and 
bk (but without the superscript "-B" ) for the Fourier coefficients of the restriction of h to the 
extremities of a cylinder: 



ike 



we study the dependence upon and bk, k & {—N, —N + 1, . . . , — 1, A^} of the measure 
on i^^, disregarding all other coefficients. The object obtained this way is a measure 
on M2(2Ar+i) concentrated on a finite set. Keeping N fixed, we then take the mesh 
a on to zero. If the limit of the measures on m2(27v+i) grists, presumably as a continuous 
distribution, the limit as the number 2N + 1 of Fourier coefficients is taken to infinity can 
be considered. We name the limiting object 

mn = lim lim m"r\^ . (7) 

This measure, if it exists, is therefore defined on a space Sji with coordinates {{ak}, {bk}) and 
we shall denote the elements of this space by 0/. This measure is to be compared with the 
probability measure induced by on the space of (p/R. (The radius of compactification 
R appears here because we normalized the jumps of the function h to be ±7r, forcing h 
to change by an integral multiple of 2it as 6 winds around one extremity.) The first, and 
principal, question is: 

(i) does the measure mn exist? 
The parallel just suggested can be pushed further. We introduce first the derivative H = 
dh/dO of the restriction of h to either extremity. It is clearly a sum of delta functions 
concentrated half way between those sites p and q at the boundary such that a{j)) ^ cr(g). 
The mass of each jump is ±7r. We shall use the letter Aj. for its Fourier coefficients. 



k& 



Clearly = ikak- At the other end we use Bk = ikbk- We will use Ak equally for the coor- 
dinates parametrizing tjjj, the derivative For the boson, the probabilistic interpretation 
of the partition function Z implies that the A;-th Fourier coefficient of the restriction is 
distributed (up to a normalizing factor) as exp(— 2A;|af p) in the limit q = 0. Consequently, 
if we use the Fourier coefficients of 

with ika^ = RCk, the probability density is e~^-^^l*^'=l^/'^, again up to normalization. For a 
long cylinder the parallel drawn here raises the following questions on m^) granted, of course, 
that the answer to (i) is positive: 

(ii) are the random variables defined by the Fourier coefficients Ak of ipi distributed nor- 
mally as e-^^l^'^l"? 




Figure 1: The distribution of JJA^, k = 1,2, 4, 8, with respect to the measure m^. The mesh 
size a corresponds to the lattices 59 x 401, 79 x 157, 157 x 1067, 199 x 397, 397 x 793 and 
793 X 1585. The curve 59 x 401 is at the top when ^Ak = 0. (See text.) 



(in) is there a constant R = Rb such that the constants f3k are simply related by 

A = ^? (8) 

(iv)is the joint distribution a product of independent single- variable distributions? 

The rest of this section will describe the response to these questions provided by numerical 
simulations. The next will provide evidence that this limit measure is both universal and 
conformally invariant. 



2.2 The distribution of h at the boundary of a long cyUnder. 

The diagrams of Figure |I] are some evidence for the existence of the distribution on the 
space Sjj. They represent the probability distribution densities of H restricted to one of the 
extremities of various cylinders, in terms of a single variable (either ^Ai, ^A2, ^A^ or ^Ag), 
all others being disregarded or, thinking in terms of the limit, integrated out. By rotational 
symmetry these densities are (almost) identical to those with respect to the imaginary part 
of the same coefficients. (A small discrepancy could arise from the fact that the numbers of 
sites along the circumference were not divisible by 4.) The square lattices contained 59 x 401, 
79 X 157, 157 X 1067, 199 x 397, 397 x 793 and 793 x 1585 sites. The first number (LV) 
is the number of sites around the circumference and is half the number of sites along the 
length (LH) minus one, or less. The Appendix gives some further technical details on the 
simulations. We note at this point that the partition function (^) is obtained by summing 
over the integer n parametrizing the linear term [a — b) = inR in both 0out and 0out- The 
analogue of this term for Ising configurations is straightforward: a configuration with exactly 
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Figure 2: The numbers Uk, k — 1, . . . ,10 ior the cyhnders, the squares and the disk. 

two clusters (of opposite signs) extending from one end to the other of the cyhnder will have 
two longitudinal jump lines. Depending on the choice of the jump across these lines, h will 
increase by 0, 27r and — 27r as 9 wraps around the boundary. Other (even) integral multiples 
of TT appear for configurations with more clusters crossing from one end to the other and the 
numbers 2mr can be used to partition the set of configurations. We have not differentiated 
the measure niD for these various classes. We should add that, for the cylinders studied 
in the present section, the configurations h whose linear term is zero are by far the most 
probable. The multiples ±2tt occurred with a probability about .0005; higher multiples we 
did not see at all. 

Even though the raw data clearly differentiate the curves attached to smaller cyhn- 
ders, smoothing helped to separate the curves between the two largest one (397 x 793 and 
793 X 1585). This smoothing was done using the kernel method with a gaussian kernel; the 
smoothing parameter was chosen according to Eq. (3.28) of [Si], in which a was taken to be 
the sample standard deviation. 

The narrowing of the gaps between the curves as the number of sites is increased is a good 
qualitative argument for the existence of the limit = liniAr-»oo liiHa^o The peaks 

of the curves go down systematically as LV and LH increase, except for the dependence 
on Ai. In this case the center of the curve for 793 x 1585 lies slightly above the center 
for 397 X 793 on the small interval (-0.05,0.05). Around = -0.05 the two curves 
cross and the curve for 793 x 1585 remains below that for 397 x 793 until approximately 
3?Ai = +0.5. Prom then on (|3?Ai| > 0.5) the two curves are so close that they cross each 
other several times, probably due to the limitation of our samples. This puzzled us and was 
checked independently by two of us. We have no explanation for it. As will be seen below 
however, the variance of the samples, a more global indicator, increases systematically over 
the spectrum of all the cyhnders considered; in particular that of 793 x 1585 is larger than 



that of 397 x 793. 
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Figure 3: The numbers Cok for 59 x 401 and 397 x 793 with the horizontal axis scaled 
proportionately to 1/ LV . 

These distribution densities are so similar to normal curves that their variances are a 
natural tool for a more qualitative assessment of the finite size effects. In order to answer 
questions (ii) and (in), we plot in Figure |^ the numbers 

r.LVxLH _ ^ 

where E^^^^^ is the square root of the variance with respect to the variable 'StA^ for the 
cylinder with LV x LH sites. (If both questions were to be answered positively, then the 
numbers k(3k = k/2Yl\ for = \im.Lv,LH^oo^k^^^^ would be a constant. Note that we 
follow the usual statistical convention of distinguishing between the theoretical value a of 
a quantity and its measured value a.) We plotted these numbers for = 1, . . . , 10 (or, 
sometimes, k = 1,2, 4, 8), together with a linear fit of these ten points for every cylinder size 
on the square lattice considered, the largest triangular and hexagonal lattices, the anisotropic 
lattice, and the 254 x 254 square and disk geometries. The latter will be discussed in in 
Section ^. The data, read from the top, appear in the order: cylinders of size 59 x 401, 
79 X 157, 157 X 1067, 199 x 397 for the square lattice Gn; of size 464 x 1069 for the hexagonal 
lattice Qhex] then of size 397 x 793 for Qo; of size 312 x 963 for the anisotropic lattice; 
of size 416 x 721 for the triangular lattice ^a; the cylinder of size 793 x 1585 for and 
the square of size 254 x 254 are superimposed; and finally the disk of radius r = 300.2 for 
Q\j. The numerical data are also recorded in Table I for k = 1,2,4,8 together with those 
for triangular and hexagonal lattices on a cylinder and those on an ellipse covered by an 
anisotropic lattice. The digit after the vertical bar gives the statistical error on the digit 
just before; for example, the first element in the table (1.609|3) means that 1/2(S^^^'^°^)^ 
is 1.609 with the 95%-confidence interval being [1.606,1.612]. The (statistical) error bars 
were not drawn on Fig. |^ as their length is approximately the size of the symbols used, or 
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less. All the linear fits meet in a very small neighborhood on the vertical axis. For the two 
cylinders with the greatest number of sites, the disk and the two squares, the ordinates at 
the origin are all in the interval [1.47071,1.47262], while the largest cylinder and the disk 
meet at essentially equal values (1.47071 and 1.47095 respectively). It is likely that, for the 
two smallest cylinders, a positive quadratic term would have improved the fit and narrowed 
the gap with the intersection of the others. 

Figure ^ reinforces this impression. The numbers tdj^^^^^ were drawn for all the linearly 
independent Fourier modes (but the constant one) for the cylinders 59 x 401 and 397 x 
793. Since each function ifssx^oi jg j-j^g g^^^ multiples of the same 59 5-functions on the 
circumference, it can be identified with a point in M^^ that we choose to parametrize with 
Aq,^Ai,'^Ai, . . . ,^A29,'^A29. Again the distributions with respect to ^A^ and '^A^ are 
identical and the corresponding samples can be united. The 29 crosses on the plot are the 
data for 59 x 401 and the 198 dots are those for 397 x 793. The horizontal axis was scaled 
differently for the two cylinders: the data were spread evenly on the interval [0, 1], starting 
at ^ for 59 x 401 and at for 397 x 793. The crosses and the dots follow almost the same 
curve when scaled that way. Hence, the change in the slopes for the various cylinders (Fig. 
0) can be seen to be the effect of calculating the slope of a curve at the origin taking 10 
values lying in an interval of length proportional to 1/ LV. This is confirmed by a log- log 
plot of these slopes (Fig. ^). The six dots can be fitted linearly and the slope is found to be 
—1.031 or, if the two smallest cylinders are discarded, —1.008. These results are indeed very 
close to —1. (It is this second fit that is drawn of the figure.) Consequently the numbers 
k/2'L\ are likely to be all equal to one and the same constant 2R\ whose four first digits are 
1.471. 

This observation together with the previous data indicates that the distribution mo 
quite probably exists and that the variances cr^ with respect to the variables ak = A^/k are 
inversely proportional to k: 

4 = = if (9) 

with est = 0.3399 close to, but unlikely to be, |. We have not discussed yet whether the 
distributions are gaussian but the form (^ is in fact in agreement with the form (|^), at least 
for the variances of the distributions with respect to one of the variables when all the others 
are integrated. The constant R\ is therefore 0.7355. 
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Geometry (lattice) 


Size 






004 


Us 




59 X 401 


1.609|3 


1.742 


3 


2.020 


4 


2.627|5 




79 X 157 


1.573 


2 


1.675 


2 


1.880 


2 


2.314 


3 


Cylinder (^n) 


157 X 1067 


1.520 


3 


1.574 


3 


1.676 


3 


1.883 


4 


199 X 397 


1.506 


6 


1.550 


6 


1.629|7 


1.793 


8 




397 X 793 


1.494|3 


1.511 


3 


1.553 


3 


1.631 


3 




793 X 1585 


1.482 3 


1.491 


3 


1.512 


3 


1.553 


3 


Cylinder (^n, Jh = 2J^) 


312 X 963 


1.487 


1.507 


1.540 


1.614 


Cylinder (^a) 


416 X 721 


1.491 


1.496 


1.536 


1.593 


Cylinder 

yr'nex J 


116 X 267 


1.599 


1.719 


1.946 


2.41^ 




235 X 535 


1.535 


1.601 


1.717 


1.952 


464 X 1069 


1.502 


1.530 


1.560 


1.716 


Disk (Qa) 


r = 300.2 


1.474|3 


1.482|3 


1.487|3 


1.506|4 


Ellipse {Qa, Jh = 2J„) 


major axis = 749.2, 
minor axis = 485.2 


1.477 


1.480 


1.489 


1.505 


Square (Qa) 


80 X 80 


1.502|4 


1.535|4 


1.600|4 


1.728 


5 


254 X 254 


1.480 5 


1.494 5 


1.510 5 


1.552 


5 



Table I: The numbers u^, /c = 1, 2, 4, 8 as measured on the cylinder, the disk, the ellipse and 
the square. Only the square lattice on cylinders is discussed in this section. See Section || 
for the others. 

We turn now to question (ii): are the Fourier coefficients Ak of ipi distributed normally 
as 6"'^*'^'=' ? To address this question we used three complementary methods that we shall 
refer to as the graphical method, the method of moments, and the method of goodness-of-ffi. 

Graphical methods seem a coarse way to assess whether an empirical distribution is a 
given theoretical one. Still they are a natural ffist choice among the arsenal of statistical 
techniques designed for this purpose. Figure |^ plots the empirical histograms for mo as 
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Figure 4: Log-Log plot of the slopes of the hnear fit of Uk as function of LV. 




-4-2 2 4 



Figure 5: Comparison of the empirical distribution ml,^^^^^^^ as a function of ^A^, k = 
1, 2, 4, 8 with gaussian densities whose variance is the sample variance. The case A; = 1 is at 
the top. 
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measured on the cylinder 793 x 1585 as functions of a single variable (]R.Ai,^A2,^A4^ and 
JRAg), all other dependence being integrated out. For these plots we have joined the data for 
^Ai,'^Ai,^Bi and '^Bi which brings the sample to 1424000 configurations. (The symme- 
tries of TTiD insure that these variables are identically distributed. We are not assuming here 
that they are statistically independent. This will be discussed in the next paragraph.) Be- 
sides these four empirical distributions, four normal curves have been plotted whose variances 
are those of the data. (These variances can be deduced from Table I.) We have left these 
empirical distributions as they are, in contrast to those seen on the Figures 0, to distinguish 
them from the (smooth) normal curves and to give to the reader an idea of the difference 
between raw and smoothed data. For the dependence on ^A2, ^A^ and JJAg, the ragged and 
smooth curves are essentially identical and, based on this evidence, one is tempted to claim 
that rriD is distributed normally with respect to these variables. The agreement between the 
two curves for ^Ai is clearly not so good. The empirical curve lies above the normal curve 
at the center, crosses it before \^Ai\ = 0.5 and remains under it at least till \^Ai\ = 1.0. 
The departure from normality is statistically significant for the dependence upon ^Ai. After 
having observed this fact, one also sees, looking more closely, a gap at the center of the curves 
for ^A2, though on a significantly smaller scale. (It might not even be visible if Figure |^ has 
been too compressed.) Since this departure from normality surprised us and, especially, as 
it is easily observable only for ^Ai, we tried to explain it as a finite-size effect. The curves 
for the smaller cylinders are however similar and the gaps seem similar to the eye. (The 
curves for ^Ai become wider as the number of sites is increased, as is seen on Fig. |I|, but the 
variance of each sample also becomes larger.) If one is convinced of conformal invariance, 
discussed in Section |^, one can also use the data from an analogous simulation performed 
on a disk whose boundary contained 2400 sites. For this geometry, the two curves for ^Ai, 
empirical and gaussian, show a similar gap. Thus, on graphical evidence only, we cannot 
conclude that the gap seen between m£i(y4i) and the normal curve is a finite size effect and 
that it is likely to disappear as LV, LH —>■ cxd. The other distributions (for A2, A^^ and A^) 
are, however, extremely close to gaussian. 

Our first attempt at a more quantitative statement is through calculation of the moments 
of the samples. We shall quickly see, however, the limitations of this approach. We denote 
by /i^j the i-th moment of the distribution with respect to the variable JJA^ 



The even moments of the normal distribution are known to be the mean (the 0-th mo- 
ment, in our case by definition), the variance cr^ (the second moment, in our case an 
unknown) and ^2s = (2s — 1)!!(T^'*. The first five non-vanishing moments are therefore 
0"^, 30"'', 15(7^, 105(7^, 945(T^°. None of the statisticians among our colleagues suggested the 
moments as a quantitative tool, probably because of the enormous errors that these mea- 
surements carry. Indeed the variance on a measurement of /x^j is {2i — l)!!cr^* if i is odd and 
((2i — 1)!! — {{i — 1)!!)^)(T^* if i is even. Consequently the error on rapidly grows out of 



hand as i increases. Nonetheless the first ten moments were calculated for the samples for 
the cylinders 59 x 401, 79 x 157, 157 x 1067, 397 x 793 and 793 x 1585. 
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Figure 6: The quotients /if j/o"! and fis^i/crl for z = 4 and 6 as functions of logL\^. 



Since the distribution is, by definition, an even function in all its variables, all the 
z-th moments, with i odd, are zero. The data only support this weakly, as more than 10% of 
the odd moments lie outside what would be the 95% confidence interval if the distributions 
were gaussian. 

Even though the errors on the moments /i^j are large, it is instructive to plot some 
of the moments as functions of logLl^, for LV = 59, 79, 157, 397, 793. Figure ^ shows the 
quotients /if j/o"! and fi^Jcr^, for i = 4,6, that should tend to 3 and 15 respectively if the limit 
distributions are gaussian. (The case z = 2 is the variance and was discussed previously.) For 
the 8-th Fourier coefficient, these quotients are monotone decreasing for both i = 4 and 6 and 
Fig. ^ repeats in another way the visual observation made from Fig. |] that the distribution 
rriD as a function of is very close to a gaussian. The plots for the first Fourier coefficient 
are less conclusive: the overall behavior is decreasing, but not systematically, and the sixth 
moment is still rather far from 15, perhaps an indication that 15 is not the limit. 

The goodness-of-fit technique is our last attempt to quantify the departure from nor- 
mality of the dependence on the Fourier coefficients, particularly of Ai. An overview of 
this technique (or more precisely this set of techniques) is given in |]dAS 
concentrate on the random variable 



We are going to 



n 



{Fn 



X] 



F{x)fdF{x) 



(10) 



known as the Cramer- von Mises statistic ( ||dAS |, chap. 4). In this expression n is the size of the 



sample and F{x) the cumulative distribution function to which the data are to be compared, 
in our case the gaussian whose variance is that of the sample. If the data Xi, X2, . . . , 



ordered (xj < Xj+i), then the empirical distribution function F„ 



X) 



Fn iS^i ■^l 1 -^2 ; 



X, 
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Figure 7: The asymptotic distribution a(w^) together with the histograms for Ai (dots) and 
As (crosses) for the cyhnder 59 x 401 and n = 1000. 

a step function defined by 

0, X < xi 
^n{x) = { ^, Xi<x<Xi+i, i = l,...,n 

The measure of integration dF{x) in w'^ is equal to f{x)dx, where f{x) is the probabihty 
distribution corresponding to F{x). The integral therefore gives more weight to intervals 
in which the random variable x is more likely to fall. This is particularly well-suited for 
our purpose as the gap between empirical and proposed distributions is precisely where the 
distribution peaks. Note that, if the data not distributed according to 

the distribution F proposed, the variable w"^ will grow with the sample size n. 

The null hypothesis Hq is, henceforth, that Fn{x) is a measurement of a variable whose 



distribution is F. Under the null hypothesis Hq, Anderson and Darling ||AD|| gave an analytic 
expression a(w^) for the asymptotic probability distribution of w^, that is, the distribution 
of the variable w'^ when the sample size n is taken to infinity. We used their formula (4.35) 
to plot the curve of Figure 0. In Chapter 4 of ||dA|j|| , Stephens indicates corrections to be 
applied to w"^ that allows finite samples to be compared to the asymptotic distribution. For 
the n's that will be used below these corrections are negligible. The two first moments of 



the distribution for are ^ (independent of n) and pS 



We concentrate on two Fourier coefficients, Ai and Ag, as our goal here is to see whether 
the departure from normality for Ai can be quantified and whether it decreases with the 
increase in the number of sites of the lattice. Again we consider 3?Aj, '^Ai, and '^Bi 
as independent and following the same statistics. We can either split the whole available 
samples into smaller sets of n elements or measure the variable w'^ for a very large n. With 
the first method, a good average w'^ can be calculated if the number N/n of smaller sets is 
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large enough. The second method will provide a single number that will, with luck, clearly 
reject Hq if it has to be rejected. We apply both. 

In splitting the large samples into smaller sets, we have to make a careful choice for n. One 
restriction comes from the actual values of the L^-integral that we want to measure. Using the 
data for the cylinder 397 x 793 in the first format described in the appendix (that is, grouped 
in 401 bins), we can estimate an order of magnitude of w^/n = {Fn{x) — F{x))'^dF{x) 
for the whole sample. (This was done using not the technique suggested by statisticians 
|dAS| , |AD|| but using rather the naive Riemann integral over these 401 bins, F[x) being 



estimated at the center of these bins. We did not attempt to evaluate the error in these 
calculations.) For Ai this integral is 0.000046 and approximately 10 times smaller for Ag. 
Even though there is an important statistical error on these numbers they give us an idea of 
the order of magnitude. We are therefore measuring a very small departure from normality 

if any. On the one hand, the strategy of splitting the sample requires to get an average 
good enough that, if it is different from |, the difference is unlikely to be of statistical origin 
and should instead indicate that Hq needs to be rejected. In other words, one should break 
the sample into several smaller samples to get a good average. On the other hand, if Hq 
is false, the quantity increases with n. Since the second moment of the distribution for 

is rather large (~ ^), we need to choose n large enough that the statistical error on 
w'^ be reasonably smaller that the number itself. A rough estimate of this error is given by 
\l A5(N/n) where N/n is the number of sets obtained by splitting the sample of size N into 
subsets of n elements. There is an obvious compromise to be struck and we chose n = 1000. 

We measured w'^ for the three cylinders 59 x 401, 157 x 1067 and 397 x 793 using the 
methods described in [|dAS| , |AD|| . These cylinders are the three runs whose data were kept 



in the second format described in the appendix, so that the exact values of all the Xi were 
available. This format allowed us to compute again the coefficients A^ (and Bi for 397 x 793). 
The two histograms for Ai (dots) and Ag (crosses) for the cylinder 59 x 401 are plotted on 
Figure ^ together with the asymptotic distribution a{w^). The range [0,0.8] accounts for 
more than 99% of the observations. Even though the crosses seem to follow more closely the 
curve than the dots, a quantitative assessment is not inappropriate. The number N/n of w'^ 
is at least 2064 for each of the three cylinders and, consequently, the statistical error on the 
resulting listed in Table II is 2^1/(45 x 2064) ~ 0.0066. Note that, for As, the intervals 
of confidence around the average w"^ always contain |, the predicted mean. Any departure 
from normality for Ag, if any, cannot be observed from this test. For Ai, the predicted | 
always falls outside of the 95%-confidence interval, though barely so for 397 x 793. This 
confirms the graphical observation made earlier and forces us to reject Hq. 

As described earlier the other strategy is to compute the numbers w"^ for a large n. We 
chose n = 250000. The disadvantage of doing so is clearly that one has a single measurement 
of w^, not an average. The results appear also in Table II. The (single) for the dependence 
on As is small for all three cylinders and the hypothesis that as the size of the cylinder (as well 
as n) goes to infinity the distribution of approaches a(w^) is totally acceptable. However, 
the values of w'^ of Ai indicate that, almost surely, they do not follow these statistics. The 
null hypothesis Hq must be rejected for Ai. 
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n = 


1000 


n = 250000 




^8 


Ai 


^8 


59 X 401 


0.1818 


0.1673 


5.578 


0.1877 


157 X 1067 


0.1802 


0.1638 


3.490 


0.0682 


397 X 793 


0.1745 


0.1656 


2.587 


0.0947 



Table II: The means w'^ for n = 1000 and the number w'^ for n = 250000. 

The null hypothesis refers, however, to a lattice of a given size and it is not these with 
which we are ultimately concerned; it is rather the limit of the distributions as the lattice 
size tends to infinity that is relevant. One obvious observation from Table II is that the 
gap between the empirical data and a gaussian curve is narrowing as the number of sites 
increases. Even though Hq has been rejected, we used the variable w^, n = 250000, to 
examine the relationship between the gap and LV. We did further runs, calculating only 
the value of for the dependence on Ai. For each of the lattices 59 x 157, 77 x 155, 101 x 
203, 125 X 251, 157 x 313, 199 x 399, 251 x 501 and 397 x 793, we obtained between 20 and 53 
measurements of the variable w^. Since Hq does not hold, we do not know the distribution 
of this random variable. On the log- log plot |], we drew the average for each lattice (x) 
together with the whole sample (dots). The spread in the sample for each lattice shows 
that the variance is very large and thus underlines the difficulty of obtaining a reliable 
mean for w"^. None the less the function id^{LV) is monotone decreasing and the linear fit 
of the data (with the first two excluded) plotted on Figure ^ indicates that a power law, 
(w^ — |) oc aLV (a ~ 2.48, e ~ —.278), is a reasonable hypothesis. We point out however 
that, with our measurements of w^, we could hardly choose between the above power law 
or any of the form {w'^ — x) cc LV^ with x in the interval [0, 1]. For this we would need 
aLV < 1/6 or LV > 10000. 

So, are the Fourier coefficients A^ distributed normally? For k large enough (say k > 4), 
it is impossible with our samples to see or calculate any difference between the empirical 
and the normal distributions. For small k, particularly for Ai, the gap is obvious but the 
goodness-of-fit technique provides clear evidence that it decreases as the size is increased. 
That the gap vanishes as LV, LH oo is not a claim on which we care to insist given only 
the present data. 



2.3 Statistical dependence and the two-point function. 

The previous paragraph studied the distribution uid with respect to a single ^Ak or '^A^, 
all the others being integrated. We now turn to the last question raised in Paragraph 
that of statistical dependence of the variables A^. and B^. 

The test for statistical independence that comes first to mind is the correlation coefficients 
between the random variables ^Ak, '^Ak, and '^Bk- These were calculated for the 

cylinder with 397 x 793 sites. According to Ch. 5 of the correlation coefficient of a pair 
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Figure 8: Log-log plot of w"^, n = 250000, as a function of LV. 
of independent gaussian variables is distributed with mean on [—1, 1] as 



0Fr( 



l-rn~dr. 11 



Here is the size of the sample used to measure the correlation coefficient (A^ = 281000 for 
the present calculation) and no longer the cutoff N used to measure the distribution m"^^ . 



If we set r = s/yN and apply Stirling's formula, (0) becomes approximately 



2n N' V2^ 

Of the correlation coefficients for a// pairs of distinct variables in SftAfc, 53^4^, '^B^ and ^B^, 
k = 1, . . . , 198, the largest turned out to be 0.0097, very small indeed. However this test 
is (almost) useless! The measure m^) is invariant under rotation of the cylinder around 
its axis, or at least, invariant under a finite subgroup. Under a rotation by an 

angle 0, the Fourier coefficient picks up a phase e**^*^ and the expected value E^A^Ai) 
must vanish unless k = —I. For pairs of variables attached to the same extremity, the 
previous numerical calculation is not useful. It is meaningful only for the pairs (JftAk, 3?-Bfc), 
{QAk,^Bk), {^Aky'^Bk) and {QAk,'^Bk) of variables at different extremities, but a more 
discriminating test of independence is certainly required. 

The two-point correlation function of spins along the boundary turns out to be a striking 
test for the independence of the variables at one end of the cylinder. Because of the identifi- 
cation cr(g) = e''^'-'^) introduced in Paragraph |2.1| , the measure m^i on the space of functions 
h, or more precisely on the space S)i, should allow for the computation of the correlation 
function {a{6i)a{62)) of spins along the extremity. Arguments have been given in the lit- 
erature (e.g. in fC^, |CZ|| ) that this two-point function should behave as the inverse of the 
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distance between the two points, namely the cord length sin((^^i — 02)/2) for the geometries 
of the disk and of the cylinder. If we distinguish between the functions h and the elements 
of the limiting space ij/, the function {a{6i)a{62)) should be (e*^'^^^!)"''^^^^))^ with 



k=l 

(The relative minus sign between the 0's removes the irrelevant constant term.) Now assume 
that the variables Ka^ and are statistically independent and normally distributed with 
variance l/(2i?sv^). Gaussian integrations lead to 



k=l 



2 



2kR% 



with 

\zk\^ = le^'^^i -e^^^^l^ = 2|l-cosfc^|, 9 = 9^ - 02. 
Since X]fcli^°s^^/^ ~ — ln(2 sinx/2), we obtain up to an (infinite) constant 

^ " sin°(^/2) 

with a = I/R'b ~ 1.360, a number that is not a simple fraction and certainly not 1. Since 
the small departure from normality discussed in the previous paragraph is unlikely to change 
much this result, there is here an obvious conflict between the prediction a = 1 and this 
result based on the hypothesis of independence of the a^'s. 

We do not know if the prediction (a (^^i) a (6^2)) oc 1/ sin((6'i —92)/2) has ever been checked 
through simulations. However the correlation can be retrieved easily from our data for the 
cylinders 59 x 401, 157 x 1067 and 397 x 793. Figure ^ presents the results together with the 
linear fits of the data after deletion of the seven first (short-distance) points. The slopes of 
these fits are 0.993, 1.001 and 0.988 for the small, middle and large cylinders. This prediction 
requires no further scrutiny. 

We are left with the possibility that the variables Ak are statistically dependent. To show 
that this is most likely the case, we offer the following two data analyses. We first study 
the conditional distributions of Fourier coefficients. Namely, we consider the distribution of 
m(JR.Ak\xmin < < Xmax), that is, the distribution of A^ when Ai is restricted to values 
between Xmin and Xmax and all the others variables are integrated. Similar conditional 
distributions with the imaginary parts are also considered. If the Fourier coefficients were 
independent, every value or interval for the restricted coefficient would lead to the same 
distribution. 

In Figure |T^, we present the distribution of ^Ai given two windows on the values of ^A2, 
for a 157 X 1067 cylinder. The windows were chosen in such a way that both distributions had 
similar statistics. The numerical data clearly show that the two distributions are different, 
and thus that these two Fourier coefficients are correlated. However, this correlation could 
be affected by the finite size of our lattices. This question of the importance of such effects 
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Figure 9: Log-Log plot of {a{6i)a{62)) as a function of sin((6'i — 62) /2) for the cylinders 
59 X 401 (top), 157 X 1067 and 397 x 793 (bottom) together with linear fits. 

is difficult to address. Since we have easy access to only three cylinder sizes, we omitted a 
rigorous study of finite-size effects. 

Nevertheless, to acquire a feel for the dependence of m{^Ak\xmin < 3^^; < Xmax) on the 
choice of /c, / and the finite size, we computed, for several values of k and /, the ratio of 
the variances of the conditional distributions of '^A^ when \^Ai\ > 1.125 and \^Ai\ < 1.125 
(which we will denote r{^Ak, ^Ai)). We also made the same comparison for the real part of 
Ak and the imaginary part of Ai. If the distributions were independent, all these ratios would 
be one. We studied these ratios for cylinders of size 59 x 401, 157 x 1067 and 397 x 793. The 
first observation is that almost all these ratios diminish when lattice size increases, so that 
there is an finite-size effect. For example, r{^Ai, '^Ai) goes from 1.19 for the 59x401 cylinder 
to 1.05 for the 397 x 793 one. Besides this finite-size effect, comparing ratios for different 
values of k and /, we observed that the statistical dependence of Fourier coefficients Ak and Ai 
diminishes rapidly when \k — l\ increases, and is weaker for larger k or /. For instance, for the 
biggest cylinder, r{^Ai,^A2) = 1.06, while ri^A^^^A^) = 1.01 and r{?R:Ai,^Au) = 0.99. 
These numbers are not conclusive, and further experiments would be essential were there 
not another more compelling argument to establish the dependence of the variables. 

As the second analysis we measure the two-point correlation (^e'^^'t'(^'^')~'t>i^^))'^ using the 
measure This is not the same as directly measuring {a{6i)a{62)) from the config- 

urations as we just did to obtain Figure Recall that rriD is obtained by the limit 
mo = limAT^oo linia^o "^D^ (see eq. (|^)). Consequently we need to set a cut-off N and 
compute the correlation function on a sufficiently large cylinder using as an approximation 
for the truncation of h to its first Fourier coefficients. If the cylinder is large enough, 
the distribution mj^ as a function of A^, k = —N, . . . ,N will be fairly well approximated by 
m°i^^ . There remains the limit N 00. To a good approximation this limit may probably 
be forgotten altogether. The previous analysis showed that the dependence between Fourier 
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Figure 10: The conditional distribution of on the 157 x 1067 cyhnder. The top graph 
contains configurations with \^A2\ > 1.125, while the lower one contains those with \^A2\ < 
1.125. 

coefficients with small indices and those with large ones is significantly smaller than the de- 
pendence amongst the first Fourier coefficients. If this is so, the gaussian approximation and 
the independence hypothesis are good ones for the distribution of A^, k > N. If the function 
h being approximated is smooth enough, the error around — 6*2 = vr for example should be 
of order o{j^) according to the computation leading to (p!2D. By definition the functions h are 
piecewise continuous and their smoothness might be improved by smearing functions as in 
the usual mathematical treatment of Green's functions. (See Section |^.) We performed the 
calculation with and without smearing. The results for the cylinder 397 x 793 are shown on 
Figure |Tl|. The thick curve is the log-log plot of {a{9i)a{92)) as a function of sin((^^i —62) /I) 
that was plotted on Figure The middle, undulating curve has been obtained by repeating 
the following two steps over the whole set of configurations: first replace the function by 
the truncation h of h to the sum of its 30 (= A^) first Fourier coefficients and, then, add 
the resulting complex number e^^K^^)-^^'^)) to the sum of the numbers previously obtained. 
Only the real part of the average is plotted as the imaginary one is essentially zero. The 
first term neglected by the truncation (031) is responsible for the wavy characteristic of the 
curve. The local extrema occurs at every 6 or 7 mesh units in agreement with the half-period 
(397/31/2 K, 6.4). A linear fit of this curve (after deletion of the seven first data) has a slope 
of —1.027. The top curve was obtained in a similar fashion, except that the two steps were 
preceded by the smearing of the function h. This smearing was done by convoluting the 
functions h with a gaussian whose variance was 2.5 in mesh units. The wavy structure is 
essentially gone. The curve appears above the two others because the smearing introduces 
in h{6i) and h{02) contributions of spins at points between Oi and O2 and thus more strongly 
correlated. The smeared correlation function is therefore larger than the two others. A linear 
fit with the deletion of the same short-distance data gives nevertheless a slope of —1.062. 
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Figure 11: Log-Log plot of ((t(^i)(t(^2)) and of {e^^m)-^))) as functions of sin((^i - 62) l^) 
for the cylinder 397 x 793. (See text.) 



The conclusion is thus that the random variables Ak (or a^) are statistically dependent 
and that the computation of (^e'^('t'^^^^~^^^^^^'j using the distribution tud leads to the predicted 
critical exponent a = 1 for the spin-spin boundary correlation. A consequence of the statis- 
tical dependence is that we cannot offer as precise a description of the measure as would 
have been possible if the answer to question (iv) had been positive. This detracts neither 
from its universality nor from its conformal invariance. 
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3 Universality and conformal invariance 
of the distributions of h on closed loops. 



3.1 Two hypotheses. 

Various crossing probabilities were measured in ||LPS|| for several percolation models at their 



critical points. Their fundamental character was stressed by two general hypotheses, one 
of universality, the other of conformal invariance, that were convincingly demonstrated by 
the simulations. The same two hypotheses will be demonstrated for the Ising model at 
criticality in Section |^. In this section, we propose similar hypotheses for the distribution of 
the function h introduced above and confront them with simulations. 

We have considered in the previous section the Ising model on the square lattice. Other 
lattices could be used. The strength of the coupling could vary from one site to another. 
Aperiodic lattices could be considered or even random ones. It is, however, easier to be 
specific and to consider two-dimensional planar periodic graphs Q. We adopt, as in | LPS |, 



the definition used by Kesten 0] in his book on percolation: (i) G should not have any loops 
(in the graph-theoretical sense), (ii) Q is periodic with respect to translations by elements of 
a lattice L in of rank two, (Hi) the number of bonds attached to a site in Q is bounded. 
(iv) all bonds of Q have bounded length and every compact set of intersects finitely 
many bonds in Q and (v) Q is connected. An Ising model is a pair {Q, J) where J is a 
positive function defined on bonds, periodic under L. The function J is to be interpreted 
as the coupling between the various sites. Only some of the models {Q, J) will be critical, 
or, as often expressed, each model is critical only for certain values of the couplings J. The 
following discussion is restricted to models at criticality. 

Let D be a connected domain of whose boundary is a regular curve and let C be a 
parametrized regular curve (without self-intersection) in the closure of D. If {Q, J) is an 
Ising model, one can measure the distribution mD c{{a'k}',Q, J) as we did in the previous 
section for m£, on the square lattice. (Although the coordinates will ultimately become 
our preferred coordinates, we continue for the moment with the a^.) The limit on the mesh 
can be taken either by dilating C and D with the dilation parameter going to infinity while Q 
fixed or by shrinking the planar lattice Q uniformly while keeping C and D fixed. As before 
we shall assume that the limit measure exists for every regular C. The previous section gave 
strong support for this supposition when C is the boundary of D and {Q, J) the isotropic 
Ising model on the square lattice. We examine the following hypothesis. 

Hypothesis of universality: For any pair of Ising models iG,J) and {G',J'), there 
exists an element g of GL{2,W) such that for all D and C 

"^D,c({afc}; J) = mgD,gc{{al}] Q\ J')- (13) 

The notation gD and gC stands for the images of D and C by g. The Fourier coefficients 
are obtained by integrating on gC with respect to 9^ , the image by the linear map g 
of the parameter 9 on C . The transformation g does not affect the underlying lattice Q. 
For example, if Q' is the regular square lattice, it remains the regular square lattice. The 
domains D and gD are simply superimposed on Q and on Q' . For the usual Ising models, 
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those defined on other symmetric graphs (the triangular and the hexagonal) with constant 
coupling or the model with anisotropic coupling on a square lattice, the matrix g is diagonal. 
It is easy to introduce models for which g would not be diagonal. We have not done so for 



the Ising models, but an example for percolation is to be found in |[LPS|| . 

To introduce the hypothesis of conformal invariance of the distributions rriD^ciWk}] Q, J), 
it is easier to restrict at first the discussion to the Ising model on the square lattice Qc\ 
with the constant coupling function Jg. A shorter notation will be used for this model: 
'>TT'D,c{Wk}) = ^D,c{{0'k}] Qn, Jo) ■ We endow with the usual complex structure, in other 
words we identify it with the space of complex numbers in the usual way. For this complex 
structure any holomorphic or antiholomorphic map defines a conformal map, at least 
locally. Given two domains D and Di we consider maps that are bijective from the closure 
of D to the closure of Di and holomorphic (or antiholomorphic) on D itself. Thus Di = (f)D. 
Let (j)C be the image of C. 

Hypothesis of conformal invariance: If (p satisfies the above conditions, then 

mD,ci{a'k}) = m^D,^ci{at}) (14) 

where the Fourier coefficients af. appearing as arguments ofm^r>,(t>c o'^e measured with respect 
to the arc-length parameter on (pC in the induced metric, or equivalently as: 



Zn 



2tt 

h"^^ o (l)c{e)e-'^^de (15) 



where (pc is the restriction of to C , h'^^ is the function h on the domain (pD and 9 is the 
(usual) arc-length parameter of the original loop C. 

Even though we have formulated this hypothesis for the Ising model on the square lattice 
with constant coupling, it is clear that it can be extended to any model {Q, J) using the 
hypothesis of universality. 



3.2 Simulations. 

Since the curve C is no longer necessarily an extremity of a cylinder, our first step is to 
acquire some intuition about the measure m^)^^^ for curves C inside the domain D. To do 
so, we continue our investigation of the cylinder for {Qui Ju)- Thus D remains the cylinder, 
but we select several curves inside it. On the cylinder 397 x 793 the curves Ci are sections 
coinciding with the leftmost column (Co), the 9-th column (Ci), the 17-th (C2), the 33- 
rd (C3), the 65-th (C4) and the middle column (C5). These curves are at a distance of 
0,0.0201,0.0403,0.0806,0.161 and 0.997 from the boundary measured as a fraction of the 
circumference. We have not checked that the measurements on curves and their mirror 
images with respect to the middle of the cylinder are statistically independent. The closest 
pair (the curves on columns 65 and 729) are, however, at a distance of 665 mesh units, that is, 
more than | of the full length of the cylinder. So to the distributions on the first five curves 
(all but the central one), we have joined those on their mirror images, doubling the numbers 
of configurations studied. Figure |T3| presents the measure rriD^Ci functions of the real part 
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Figure 12: Two "typical" configurations on a disk of radius 200 witli free boundary. 




Figure 13: Distributions m£)^Ci as functions of the real part of in the natural order: = 1, 
k = 2,k = A,k = 8,k = 16,'/c"= 32. 



of the Fourier coefficients A^, k = 1,2,4,8, 16,32. Each graph shows the dependence on a 
fixed Ak for the six curves. On each graph the lowest curve at the origin corresponds to Co = 
boundary, the case studied in Section ^ As the curve C is taken closer to the center of the 
cylinder, the distribution becomes sharper at the center. This is perhaps to be expected as 
the sites at the boundary are freer to create clusters of intermediate size than are the sites 
in the bulk, increasing thereby the values of the various Fourier coefficients. (See Figures 



12 and IS.) Another natural feature is the gathering of the distributions for all the interior 



curves on the plot for ^Aiq and 3^^432. Indeed the higher Fourier coefficients Ak probe small 
scale structure, at the approximate scale of ^ in circumference units. For example, the 
Fourier coefficient A32 will be sensitive mostly to clusters having a "diameter" of ~ 12 mesh 
units or less and these clusters intersecting the curves C at a distance of 32 or of 64 mesh 
units from the extremity should be distributed more or less the same way. In other words 
the bulk behavior is reached closer to the boundary for higher Fourier coefficients. One last 
observation about these plots is that the bulk distribution is definitely not a gaussian in 
^Ai\ It is sharply peaked at the center but still has a wide tail. (The distribution in ^Ai 
measured along the mid-curve of the cylinder can be better seen on Figure |14| below.) 



27 



To examine the hypotheses of universahty and conformal invariance, we ran simulations 
on other pairs {Q, J) and on other geometries. We discuss both at the same time. Three 
other pairs {Q, J) were considered: the regular triangular and hexagonal lattices Q/\ and ^hex 
with the constant function J and the regular square lattice with a function J that takes 
a constant value Jh on the horizontal bonds and another constant value on the vertical 
ones with = 2Jy. We shall call this model the anisotropic Ising model. This choice of 
J makes the horizontal bonds stronger than the vertical ones and clusters of identical spins 
will have a shape elongated in the horizontal direction as compared to those of the isotropic 
model (^n, Jn). 

The critical couplings are determined by sinh 27^ sinh 2 J„ = 1 (see |B[ or |[MW|| ). If the 
hypothesis of universality is accepted then it follows from formula (5.9) of XI.5 of ||MW|| that 



the matrix g that appears in (|13]) (when (Q, J) is the critical model on the square lattice 
and {Q', J') the anisotropic model) is0 

1 
sinh 2Jh 

(The critical value of Jh for which Jh = 2J^ is 0.609378....) The lattice used for the 
anisotropic model has LV = 312 and LH = 963. These dimensions correspond to a cylinder 
on the square lattice with a horizontal/vertical ratio of 1.999, very close to the one used 
for the square lattice 397 x 793 that has a ratio 1.997. The lattice used for the triangular 
lattice was oriented in such a way that every triangle had one side along the horizontal axis 
and the dimensions used were LV = 416, that is, the number of horizontal lines containing 
sites, and LH = 712, the number of sites on these lines. The aspect ratio for a square 
lattice corresponding to these numbers is 2.001. The largest hexagonal lattice used was of 
size 464 x 1069. Again 464 is the number of horizontal lines containing sites and LH = 1069 
is the length of the cylinder in mesh units. The corresponding aspect ratio for the square 
lattice is 1.995. We also measured the smaller hexagonal lattices of sizes 116 x 267 and 
235 X 535. The difference between these four ratios is smaller than the limitation due to 
finiteness discussed in [[LPP|j|]. The distances of the curves Cj from the boundary were chosen 



as close as possible to those used for the cylinder on {Qa, Ja) and given above. (The manner 
in which the Fourier coefficients of the restriction of h to these curves were calculated is 
described in the appendix.) 

As evidence for the hypothesis of conformal invariance, we compared three different 
geometries, namely the cylinder used in Paragraph a disk, and a square. We identify the 
cylinder with the rectangle in the complex plane of height v (its circumference) and of length 
h. The analytic function z ^-"^-^z/v j^^g^pg ^^j^jg cylinder onto an annulus. With our choice 
of dimensions for the cylinder {v = 397, h = 793), the ratio of the inner and outer radii is 
less than 10~^ and unless the outer diameter of the annulus is larger than 10^, the inner 
circle contains a single site. We took the liberty of adding this site to the domain and of 



^The point is that because of the anisotropy the two-point correlation function decays more slowly in the 
horizontal direction than in the vertical, behaving at large distance as l/(x^ + a'^y'^y/'^ with a = sinh2J;i « 
1.54. The appropriate conformal structure is that defined by the ellipse + a^y^ — 1. We are grateful to 



Christian Mercat for this reference. See also his thesis (|Me|) in which the conformal properties of the Ising 
model are discussed from quite a different standpoint. 
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identifying it with a disk. In other words, although the geometries of the cyhnder and of the 
disk are not conformally equivalent in the sense of the hypothesis, the finite size realization 
used here for the disk differs by a single site from the annulus conformally equivalent to the 
cylinder. The radius of the disk was taken to be 300.2. The disk can be mapped onto the 
square by the Schwarz-Christoffel formula 



which defines a map, with the unit disk as domain, holomorphic except in the four points 
_l_g±i7r/4_ gQi;]-^ ixiaps satisfy our requirements. For the square and the disk, the distributions 
were measured at the boundary. For the disk, they were also measured on the four circles 
corresponding to the inner circles on the cylinder that are not at its center. This latter 
circle on the cylinder is mapped, inside the disk of radius r = 300.2, onto a circle of radius 
~ 0.57, less than one mesh unit. The distribution mo,c on this circle is clearly impossible 
to measure for this lattice size. 

Table I of the previous section has been completed with the data ook = k/2{llkY for 
six new experiments: the three new Ising pairs {Q, J) (triangular lattice, hexagonal lattice 
and anisotropic function J) on the original cylindrical geometry; the two new geometries 
(disk and square) covered by the square lattice; an ellipse covered by the square lattice with 
anisotropic interaction. For the square, two runs were made on a lattice of 80 x 80 and 
254 X 254 sites. The data for the disk and both squares were also drawn on Figure |^. As 
discussed previously, it can be seen there that their cj^'s follow exactly the same pattern as 
those of the cylinders and that the ordinate at the origin of their fits falls in the same very 
small window [1.47071,1.47262]. It is interesting to notice that the small lattice 80 x 80 
on the square geometry leads to tu/j's that are between those of the lattices 199 x 397 and 
397 X 793 for the cylinder. Considering that the number of sites in the lattice 199 x 397 is 
more than twelve-fold that in the small square, this might seem surprising. The explanation 
is likely to be that the number of sites on the boundary where the distribution is measured 
is the leading cause of the finite size effect. The cj^'s for the triangular lattice and for the 
anisotropic model were obtained from the 401-bin histograms of the empirical distributions. 
(See the appendix.) No attempt was made to provide confidence intervals. The linear 
fits of the ojk.k = 1,2,4,8, are 1.4723 + 0.0152A; (triangular lattice) and 1.4695 + 0.0180A; 
(anisotropic model). For the largest of the square lattices it was 1.4712 + 0.0102/;; and for the 
disk 1.4710 + 0.0044A;. The ordinates at the origin (1.4723, 1.4695 and 1.4712) are extremely 
close to the narrow window above for the larger cylinders, the disk and squares, especially 
striking as the samples for these experiments (200K) were the smallest of all in this section 
and the previous one. The linear fit for the largest of the hexagonal lattices is 1. 4793+0. 0294A; 
and the ordinate at the origin is not quite so good but the slope remains large compared 
with the other fits. Indeed the product of the slope and the circumference LV is in the 
four cases: ~ 8.1 (square); ~ 6.3 (anisotropic); ~ 6.4 (triangular); ~ 14. (hexagonal). This 
suggests that the circumference of the hexagonal lattice must be twice that of the triangular 
lattice in order to obtain comparable results, perhaps because it contains only half as many 
bonds per site as the triangular lattice. The ordinate at the origin (1.4770) is nevertheless 
close and this is important because it confirms the suitability of the construction of the 




dw 



(16) 
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Figure 14: Distributions mD,c'i functions of the real part of Ak {k = 1 on the first hne, 
k = 8 on the second) on three different curves C: the boundary (first column), the curve 
at a distance of 16 mesh units on the cylinder 397 x 793 and its conformal images (second 
column) and the curve at the middle of the cylinder (last column). 

function h that is described in the appendix, a construction less obvious and more difficult 
to implement for the hexagonal lattice than for the others. The anisotropic lattice on an 
ellipse was included to demonstrate that the measure m^.c at the boundary is able to select 
the appropriate conformal structure even when it is not obvious by symmetry. One map 
between the structure attached to the anisotropic lattice, thus the square lattice with the 
indicated asymmetric interaction, and that attached to the square lattice with symmetric 
interaction takes an ellipse x"^ + ay"^ < 1, a = 1.54369, to the disk + < 1. As ellipse we 
took one whose major and minor axes were of lengths 749.2 and 485.2. The usual linear fit 
of ojfc, = 1, . . . , 10 yielded 1.4712 + .0044x. 

The plots of Figure [1^ show the measure as a function of and ^Ag, when h is 
restricted to three different curves on the cylinder or to their conformal images on other 
geometries: the boundary, the second inner circle (at a distance 0.0403 from the boundary 
measured as a fraction of the circumference) and the circle in the middle of the cylinder. For 
the boundary (first column of Figure |1^ five models have been drawn: the cylinder covered 
with the square, the triangular and the anisotropic lattices, the disk and the square both 
covered with the square lattice. (The numbers of sites on the various lattices are those given 
earlier in this paragraph; only the data for the square of side 254 were drawn here.) For the 
second column of the figure, the same models were used but no measurements were made 
for the square. For the curve in the middle (third column), only the three lattices on the 
cylinder were measured, because the corresponding circle on the disk is too small to allow 
for reliable measurements. (See below.) To these three lattices a fourth square lattice, with 
199 X 2399 sites, was added on the ^Ai plot. The agreement is convincing, as it is for the 
distributions along the other curves Ci that we measured. 

At first glance no cogent comparison can be made between the central circle on the 
cylinder and a circle in the disk. A circle in the middle of a short cylinder is equivalent to a 
circle in an annulus, but when the cylinder becomes extremely long, it is more like a circle 
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Figure 15: The cufc, = 1, . . . , 8 on the boundary and on three inner circles of the disks of 
radius 100.2 (+), 200.2 (x) and 300.2 (•). The continuous hues represent the corresponding 
data for the cyhnder 397 x 793. 

in the plane. For example, if the cylinders of size 397 x 793 and 199 x 2399 are mapped to 
an annulus of outer radius 1, the inner radii will be 4 x 10^^ and 1 x 10^^^ and the images 
of the central circles will have radii 2 x 10~^ and 4 x 10~^^ respectively (too small to make 
a measurement). All circles in the plane are, however, conformally equivalent. So we can 
still compare the distributions on the central circle of a cylinder with the distribution on a 
circle in the plane. This is easier said than done, because the larger the circle the larger the 
domain needed to make useful measurements. There is none the less a method, so that the 
distributions measured on the central circle on both cylinders 397 x 793 and 199 x 2399 can 
be considered as distributions in the bulk. 

We first calculate the u^s for progressively smaller circles inside disks and observe that 
they do tend toward a limiting distribution. The main difficulty is again the finite-size 
effect revealed in Figure ^ We compare corresponding inner circles on disks of radius 
100.2, 200.2, 300.2. On each of these, the distributions were measured on inner circles of 
radius 1., 0.5, 0.4, 0.3, 0.2, 0.1319, 0.1 and 0.04790 times the outer radius. The smallest inner 
circle on the disk of radius 100.2 has a radius 4.8 in mesh units. Finite-size effects will be 
indeed important! Though we measured the uj^s for k up to 32, the overall behavior is clear 
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Figure 16: The cu^, k E {1,2,..., 32} on a disk (+) and on a cylinder (• and ■) for various 
inner circles. 

for /c = 1, . . . , 8, as presented on Figure |1^. Only the circles of relative radius 1., 0.4, 0.2 
and 0.0479 were retained for ease of reading. The "+" are for the disk of radius 100.2, the 
"x" for 200.2, the "•" for 300.2 and the corresponding data for the cylinder 397 x 793 are 
joined by straight lines. For the boundary, the three disks give a better approximation of 
the limiting distributions than the cylinder but for the inner circles the roles are exchanged. 
The spread between the three disks, and between them and the cylinder, is particularly 
important for the smallest inner circles but the way it decreases with the increase of the disk 
radius supports the hypothesis that a common distribution for these two geometries exist 
on each of these circles. 

A comparison of the distributions on inner circles for the disk and the cylinder is therefore 



possible. Figure |T6] shows the w^'s, k G {1,2, .. . ,32}, for the cylinder 397 x 793 (•) and 
the disk of radius 300.2 (+). Curves were added to help the eye. Seven circles were used. 
Their distance from the boundary of the cylinder, in mesh units, and their relative radius for 
the disk (in parenthesis) are 0(1.), 8(0.881), 16(0.776), 32(0.602), 64(0.363), 128(0.132) and 
192(0.0479). The measurements on the central circle of the cylinder were added (■). Only 
for the two smallest circles of radius 0.132 and 0.0479 is the agreement less convincing but, 
again, the previous figure showed how the gap diminishes as the outer radius increases. We 
shall therefore refer to the limiting distribution, approached by that on central circles of 
cylinders and on very small inner circles of disks, as the bulk behavior irrespective of the 
global geometry. 

Another way to check that the bulk behavior is almost reached in the middle of the cylin- 
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der is to compute the spin- spin correlation along the central parallel. Using the conformal 
map from the (infinite) cylinder to a disk one can see that the function {a{9i)a{92)) should 
be proportional to sin~°''""^((6'i — 62) /2). The conformal exponent Obuik is | (see [ [MW|| ). A 
log-log fit of (cr(6'i)cr(6'2)) as a function of sin((6'i — 62)12 gives a slope of —0.257. We can also 
verify that the measure m£, on the distributions on this curve allows us to recover this Obuik 
by the measurement of (^e^W^^')~^(^^^^'^, As in Paragraph |2.3| we did this by first smearing the 
functions h with a gaussian and then truncating their Fourier expansion at N. (As before 
we set = 30 and the variance of the gaussian to 2.5 mesh units.) The linear fit of the 
log-log plot leads to an abuik ~ —0.260, in fair agreement with the expected value. 



Although the random variable 'SiAi is not normal, Figure |1J and, less clearly. Figure [13 
show that higher Fourier coefficients are close to normal. In fact, starting around k = 8, 
the histograms of the are graphically undistinguishable from the normal curves whose 
variances are those of the samples. One may ask quite naturally if the distribution of these 
variables is given, at least asymptotically, by the law (||) with, maybe, another constant -Rbuik- 
The linear fit for the wj?""^ is 4.380 + 0.0065fc and the slope is slightly smaller than that of 
the boundary, an indication that finite-size effects are smaller in the bulk. It therefore 
seems likely that the variables are asymptotically normally distributed as in (|^) with 
(i?buik)^ = 2.190. The ratio (i?buik)V^| is 2.98, very close to 3. 

The existence of a nontrivial bulk behavior on curves in the plane, was by no means 
initially evident and may, in the long run, be one of the more mathematically significant 
facts revealed by our experiments. One supposes that the distribution of spins in a fixed, 
bounded region for the Ising model on the complete planar lattice at criticality on a lattice 
whose mesh is going to is such that they are overwhelmingly of one sign with substantially 
smaller islands of opposite signs and that these islands in their turn are dotted with lakes 
and so on. This is confirmed by the two typical states of Figure |T2| in which the large islands 
of opposite spin appear only in regions infiuenced by the boundary. Typical states for the 
cylinder are similar (Figure ^) but the conformal geometry is such that the bulk state is 
reached closer to the boundary. The conclusion is not, apparently, that in an enormous disk, 
thus in the plane, the integral of h against any fixed smooth function on a fixed smooth 
curve is generally very close to 0, so that the distribution of each of the Fourier coefficients 
^Ak and QAk approaches a (5-function. Rather, they are approaching a distribution which 
is not trivial but is, at least for k small, clearly not a gaussian. What we may be seeing is 
the effect of the shifting boundaries of the large regions of constant spin. Once a circle in 
the plane is fixed, the boundary between even two very large regions of different spin can, 
as the configuration is varied, cut it into intervals of quite different size. 

Indeed the existence of a nontrivial limiting measure on the space of distributions on the 
boundary was itself not certain beforehand. In spite of the attention we gave in Section ^ to 
the possibility of its being gaussian for the boundary of a circle, the exact form is perhaps 
of less mathematical significance than its universality and conformal invariance. 



3.3 Clarification. 



In order not to encumber the initial discussion with unnecessary abstraction, we worked 
with the distributions Tnj:, c{{0'k}, G, J)- A better theoretical formulation would be in terms 
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of a measure m^i ,7 on the set of real- valued distributions in the sense of Schwartz on the 
oriented smooth curve C, or if C were merely regular (thus sufficiently differentiable) on 
some Sobolev space. To be more precise, the measure is on the set of distributions that 
annihilate the constant functions. (To be even more precise, this is so only if the curve is 
contractible. For other curves, such as the circumference of a cylinder, the set of distributions 
whose value on the constant function 1 lies in {2rmT | m G Z, m 7^ 0}, may have a nonzero 
measure. Under many circumstances, it is small enough that it can for numerical purposes be 
supposed 0.) To introduce the measure mD,c concretely, we need a basis for the dual space, 
thus in principle just a basis for the smooth (or regular) functions on C modulo constants. 
If this basis is {(pk\k = 1, 00} then A {^{fk)} defines a map of the distributions into a 
sequence space and a measure is just a measure on the collection of real infinite sequences, 
{fii, /i2, . . . }. It would have to be defined by some sort of limiting process from measures on 
M^. The simplest such measures are product measures. Given such a measure on the space 
of real sequences, it defines, at least intuitively, a measure on distributions if for almost all 
sequences {/Xfc}, the assignment ipk — > /ifc extends to a distribution, thus in particular if it 
lies in some Sobolev space. For example, if a parametrization x(^^), < < 27r of the curve 
has been fixed then one possible choice of the basis is the collection 

{xi9) ^e'^'^xie) '^e'^^'^lk > 0} 

Then it is better to put the /i^ together in pairs and to use sequences {A^} (or a^. = A^/ik) 
of complex numbers. This has been the point of view of this paragraph. Starting with a 
given parametrization, we examined the joint distributions of the complex random variables 

The parametrization also allows us to introduce the measure dO/^ir and thus to identify 
functions with distributions. In particular, in order for a measure on sequences to yield a 
measure on distributions it is necessary, and presumably usually sufficient, that the sum 

00 

k=—oo 

converges as a distribution for almost all sequences a^. For example, if the measure is a 
gaussian defined by 

00 

exp(-a^/c|afcp) (17) 

fc=i 

then the expectation E(|afcp) is l/{2ak), so that converges almost everywhere. 

As a result, the sum (0) converges almost everywhere as a distribution. This conclusion 
remains valid provided only that the expectations {a^ai) are those of the gaussian (0), a 
property that according to the results of Section ^ the measure ttid has a good chance of 
possessing. Therefore, if A = J2k '^kO-k is any random variable that is a linear function of the 
Ofc then the expectation E{\^) is calculated as though the measure were gaussian. 

Our method is numerical, so that we approximate the measure on sequences from a large, 
finite scattering of functions /i, or rather of their derivatives, because the derivative H is well 
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defined as a distribution, although h itself is not. The distribution if is a sum of 5-functions, 
with mass ivr at each point where the curve crosses a contour line of the function H. Since 
its value on the constant function 1 is the sum of those jumps, this value is 0, as noted, 
whenever that sum necessarily vanishes, either because the curve is contractible or because 
the cylinder is extremely long. 

Although our construction required a specific parametrization, the resulting measure on 
distributions may be independent of the parametrization and, more generally, even of the 
choice of basis. We did not attempt to verify this. It may be useful, however, to describe an 
example. 

When D is a disk of radius 1 with the boundary C parametrized in the usual way by 
arc length 9, the function h can be recovered by integration with respect to d9 from the 
distribution H. The measure on distributions on C is, as we discovered, not equal to the 
gaussian measure associated to a constant, 2R\, times the Dirichlet form Q.{H), but, if we 
ignore the reservation expressed at the end of Paragraph |2.1| , the variance of linear functions 
of the Fourier coefficients can be calculated as though it were. We recall that to calculate 
or £l{H,H) if we want to stress that it is a quadratic form, we extend the function 
/i as a harmonic function to the interior and then 
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£1{H) =^{h)=^{h,h) = -J {{^r + i^ndxdy. 



or, extending it to an hermitian form. 
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if we use again the symbol h for the harmonic function inside D. 
distributions with functions by means of the bilinear form 



If we identify formally 
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2^ 
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hi{9)h2{9)d9 = {hi,h2), 



or, in complex terms. 



^ r hi{9)h2{9)d9 

Jo 



(^1,^2), 



(19) 



and regard therefore £3 and S) as operators, so that il{H) = {0.H, H) and T){h) = {'Dh, h), 
then, as a simple calculation with the functions e^^^ shows, 2) has as an eigenvalue of multi- 
plicity one, eigenvalues ^, |, |, . . . , each with multiplicity two, Q. has eigenvalues |, |, . . . , 
each with multiplicity two and on the domain of O, the orthogonal complement of the con- 
stant functions, 4D = £l~^. More precisely, and this is the best form for our purposes, if the 
Fourier expansion of h is a^e^^^ then 
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or if h is real, 



fc>0 

Suppose now that D' is any domain, C its boundary, and ip' any smooth function on C. 
The function (p' defines a hnear form 

A ^ X{ip') (20) 

on distributions. By conformal invariance, the measure on distributions on C is obtained by 
transport of the measure on distributions on C using any conformal transformation from 
D to D'. If the measure on the distributions is in fact well-defined, independently of any 
choice of basis, then the characteristic function of (^) is formally calculated as 

J exp (-2/?|(nA, A) + ia\{ip)) / j exp (-2i?|(0A, A)) 

which is 

exp(-a'n(A<^)/8i?|) 

ii ip = if' o (j) and A<^ is the distribution such that A((y9) = n(A, A^^). Consequently, the 
probability distribution of the random variable defined by (pOf ) will be gaussian with variance 
T? given by 

1/2T? = 2Rl/n{X^). 

But = il~^ip so that 

1/2S2 = i?|/2D((^). (21) 

Suppose for example that D' is a square of side | and that we parametrize its boundary 
by arc length: s = s{t), s{0) being one of the vertices. We write t = t{s) for the inverse 
function. The Schwarz- Christ off el map of the disk onto the square is depicted in Figure IT? 



where the curves intersecting at the center of the square are the image of the rays on the disk. 
Although arcs of equal length on the circular boundary are mapped to intervals of different 
lengths on the edge of the square, this effect is important only close to the vertices. Thus 
if (f'f. is the function s(t) cos{kt) and ipk = v^'^ o then the distribution of the random 
variable defined by (pk should be gaussian with variance D{ipk) / R%- Moreover 2)((/9fc) is 
obtained from the Fourier coefficients of (p>k and they are calculated by observing that, apart 
from a constant factor, which is unimportant, the Schwarz-Christoffel transformation ( [T6| ) 
restricted to the boundary w = e*^ of the unit disk can be expressed in terms of the elliptic 
integral of the first kind, 

dip 



Fm 



'1-2 sin^ iIj 



With our choice of s(0), the function s(t), for t G [0, |], is 
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Figure 17: The Schwarz-Christoffel map from the disk to the square. The images of the rays 
of the disk are the curves intersecting at the center of the square. 



The graph of (t, s{t)) on [— tt, tt] is obtained from that on [0, ^] by translation by (|, f)- The 
function s = s{t) is odd and composition of odd or even functions with s preserves their 
parity. If the basis {s{t) — > coslt,s{t) smlt,l > 0} is chosen and (pi{t) — cos{ls{t)) (or 
sm{ls{t))) written as 



For a given / the random variables ^Ai and '^Ai are identically distributed on the disk, at 
least when the number of sites at the boundary is a multiple of 4. On the square they were 
shown to be also identically distributed, at least in the limit of the simulations, when these 
variables are measured with respect to the induced parameter. However, if the arc-length 
parameter s is used on the square, the two variables ^Af and '^Af arc identically distributed 
only when / is odd. The graphs of two functions cos(/0(t)) and sin(/0(t)) are translations of 
each other when / is odd but not when / is even. The variances of the random variables ^A^ 
and ^Af must then be distinguished and they are given by 




(22) 



k>l 



then the Dirichlet form is 




(23) 



k>l 




and 




(24) 



Using these formulas we shall compute the numbers 





and o;^ 
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introduced in Section ^. 
The coefficients 



Cik = — / cos I s{t) cos kt dt and Sik = — / sin I s{t) sin kt dt 



are therefore needed. They can be calculated numerically. The convergence rate of 
is however slow. The elliptic integral F{t\2) = Jq{1 — 2 sin^ tjj)~^dilj behaves like {t — j)^ 
a.s t —>■ j~ . Consequently, the function s has a similar behavior at integer multiples of | 

3 

and the absolute values \Cik\ and \Sik\ decrease approximately as M/fca. Even with the 250 
ffist Fourier coefficients Cik, k = 1, . . . , 250, the Parseval identity for the function cos/0(t) 
is satisfied to only five decimal digits. We decided nonetheless to restrict the sums to 
these 250 ffist coefficients. Since all the terms in 2? are positive, the truncated sums will 
lead to larger estimates of the cu's than the true sums. 

Since we wish to compare cu^^^ and w^^^ with those measured with the simulations done 
on the square of side 254, it is appropriate to modify slightly the Dirichlet form ( ^31) to 
take into account finite-size effects. Paragraph |2]^ showed that the quantity /c/(2E^) is not 
strictly constant on a finite lattice but grows slowly. We found that 

k 



2Rl{l + ek) 



was a good approximation. (See Figure 
the variance, we decided to replace it by 



Since the ratio k/AR"^ in plays the role of 



4i?|(l + efc)' 



The slope e is that of the linear ffi appearing in Figure ^ for the square with 254 x 254 sites. 
Table III lists the values of uf, that is w^^, and w^^^, / = 1, 2, 3, 4, 5, as measured by the 



simulations and the values tu^^^ and uj^j^^ obtained using the (truncated) sums ( p4D with 
finite-size effects introduced as discussed. The original values ui have been added to give an 
idea of the discrepancy that the use of the arc-length parameter introduces. The values il^f 
and uf are close to one another and the latter are always greater than the former, probably 
because of the truncation. 



/ 








1 


1.480 


1.380 


1.388 






2.241 


2.251 


2 


1.494 


0.963 


0.974 


3 


1.505 


1.365 


1.386 






1.921 


1.958 


4 


1.510 


1.146 


1.181 


5 


1.520 


1.451 


1.499 



Table III: The numbers a);, oof and ujf for I = 1,2,3,4,5 for the square of side 254. When 
they differ, the u for ^Ai is placed above the uj for QAi. 
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3.4 Conditional probabilities. 

Suppose that the curve C of the previous section is the disjoint union of two curves Ci 
and C2. Then the space of distributions on C is the product of the space on Ci and the 
space on C2. We fix the model to be the Ising model at criticality on the square lattice and 
denote a distribution by ij) and the measure whose meaning was clarified in the previous 
paragraph by mo^cW- Then, in principle, the conditional probability mD,c{4'i\'4'2) on the 
set of distributions on Ci is defined for each distribution -02 on C2. Whether this is so is not 
so easy to test experimentally. To approximate the conditional probability numerically with 
our methods we have to choose a neighborhood U of 1^2 and proceed as before, eliminating 
from the sample all distributions ip' = {1^1,1^2) for which ip2 does not fall in U. We recall 
that is a distribution given by a sum of (5-functions on Ci. First of all, the neighborhood 
f/ is a neighborhood in an infinite-dimensional space, so that it is going to be, in any case, 
very large. Secondly, we cannot eliminate too many distributions for then the samples would 
be far too small. Thus U is going to have to be enormous. The notion seems nevertheless 
to be workable even at a coarse experimental level. 

There are two properties that one might expect. We can introduce and study experimen- 
tally the measure on the distributions on Ci obtained when the spins on C2, or in a small 
neighborhood of it, are all taken to be +1. This of course presupposes some kind of compat- 
ibility of C2 with the lattice structure, as in the examples studied where C2 passes through a 
row of sites, or some way, either theoretical or practical, of specifying the neighborhood, but 
granted this, we consider the measure ^«d,Ci (^^11^2; +) obtained from this familiar condition. 
It is defined quite differently than the conditional probability "^D,c(V^l|0) for ip2 = 0. (See 
Paragraph [4.2| .) None the less, one could hope that they were equal. The experiments to 
be described are too coarse to establish this with any degree of certainty, but do render the 
expectation plausible. 

The second property is the markovian property. Suppose that Ci is the disjoint union of 
C3 and C4, so that ipi is a pair {ip^, ip^). Suppose moreover - this is the essential condition - 
that C2 separates C3 from C4. Then one can hope that conditioning the measure w^d,c(V'i 1^2) 
on ip4 leads to a measure mj:i^c{4'3\4'4\4'2) that is equal to mD,c{4'3\i'2), thus the measure on 
the distributions on C3 when the distributions on C2 and C4 are given is independent of the 
distribution on C4. The influence of the distribution ip^ is not propagated across C2 when 
the distribution on C2 is fixed. 



We begin by examining Figure [T^ in which three typical states are shown, from left 
to right: free boundary conditions on a cylinder of circumference 199 and length 31; free 
boundary conditions on a cylinder of circumference 199 and length 399; and free boundary 
conditions on the left but + boundary conditions on the right on a cylinder of circumference 
199 and length 31. The picture in the center is the familiar one: towards the middle there is 
a tendency to form very large clusters of constant sign, indeed there is only one very large 
(white) cluster but at the boundary the clusters are smaller. Recall as well that for a cylinder 



there is conformal distortion. In Figure 12 the phenomenon is illustrated without distortion: 



there is one large (white) cluster on the left and one large (black) one on the right. In the 



picture on the left of Figure 18, the freedom to form smaller clusters is reinforced by the 



proximity of the two boundaries. There is almost no bulk behavior at all. On the other 
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Figure 18: Configurations on cylinders. 



hand, in the picture on the right, the boundary condition is forcing a single large cluster on 
the right and this cluster is attempting to envelop the left boundary as well. 

This qualitative description is confirmed by a calculation, for the measures on the distri- 
butions on the left boundaries, of the numbers Uk introduced in Paragraph p.2|. The results 



are plotted in the diagrams of Figure ^ for the measure associated to the left boundary in 
the three cases. In clockwise order from the upper left, they are: free boundary conditions 
on a cylinder of size 199 x 31; free boundary conditions on a cylinder of size 199 x 399; 
boundary conditions on the left free, those on the right constant, and size 199 x 31. In 
the diagram on the lower left, they are superposed. The graph in the upper right is like 
those of Figure ^, except that we have used new statistics with a smaller sample, so that 
the graph is somewhat irregular. All graphs are pretty much the same except for the first 
four or five values oi k. As far as the higher values of k are concerned the two boundaries 
are effectively at an infinite distance from each other. For k = 1, there is a pronounced 
difference between the graphs so that the distribution of on the short cylinder is flatter 
than on the long cylinder. On the other hand, when the boundary condition is imposed the 
value of u)fc increases and the distribution of ^Ai is peaked. The superposition of the three 
curves is shown in Figure 

We now take Ci to be the left boundary of a cylinder of aspect ratio 199/31 and C2 to be 
the right boundary. To test the assertion that m£),Ci('^i|C2, +) is the conditional probability 
m£) ci^pilO), we thermalize for free boundary conditions at both ends of a cylinder of size 
199 X 31 but only keep those samples for which 

< .2^2, l^Asl < .35^3, (25) 

< .2V2, I'^Asl < .35^3. (26) 

About 3 out of every 10,000 states satisfy this condition. So our crude experiments will 



\^Ai\ 


< .125, 


l^A^l 


\QAi\ 


< .125, 





not permit a substantially smaller neighborhood of 0. In Figure ^ we plot the resulting 
collection of uj^ together with those obtained from the previous experiment with + boundary 
conditions on the right side. We see that in spite of the large size of the neighborhood, the 
two graphs are quite close. It is the values of ui and uj2 that tell. The graphs of the 
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Figure 19: The numbers Cj^ associated with the conditional distributions. (See text.) 




Figure 20: The distribution of "^Ai for the three pairs (domain, boundary conditions) de- 
scribed in the text. At the center the curves are in the order, from top to bottom: short 
cyhnder with constant spins on the right, long cylinder, short cylinder with both sides free. 
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Figure 21: The numbers Uk for mD,Cii4'i\C2, +) and mD,ci4'i\0)- 
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Figure 22: The distribution of for mD^Ci{4'i\C2, +) and mc^cii^ilO)- 

distributions of are compared in Figure ^ to ensure that not only are the variances 
cfose but also the probability measures themselves. Without being at all conclusive, the 
experiment encourages the belief that 

mD,c{i'i\0) = mD,Ci{i'i\C2,+)- 

In order to test whether the probabilities are markovian we considered on the one hand 
a cylinder of size 199 x 31 on which we thermalized, keeping only the distributions that on 
the right boundary satisfied the conditions 

\^Ai\ <.125, \^A2\ < .2V2, l^A-il < .35^3, 

(27) 

• 3 < QAi < 1, \QA2\ < .2V2, I'^A-il < .35^3. 

On the other hand we considered a cylinder of size 199 x 61 on which we thermalized with 
spin + as the boundary condition on the right and then selected only those states satisfying 
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Figure 23: Test of the markovian hypothesis. The numbers Cjk for the two cyhnders with 
3?Ajt on the left and '^Ak on the right. 




Figure 24: The distributions on the left boundary as a function of 3?Ai (on the left) and of 
(on the right). 

the conditions ( ^Tf ) on the distributions for the central meridian. We then examined the 



resulting measure on the distributions on the left boundary, in particular the distribution of 
9fJy4i and ^Ai. The markovian hypothesis asserts that, when we fix the distribution on the 
center, the measure on the distributions on the left boundary is completely shielded from the 
boundary conditions on the right, although once again we are prevented by the necessity of 
allowing the rather large open neighborhood (^) from actually fixing the distribution on the 
center. We can only impose very crude constraints on the first few Fourier coefficients. For 
the experiments on the smaller cylinder about 3 samples in 10,000 are kept; on the larger, 
curiously enough, about 1 in 1,000. 

In Figure |2^ the Cjk are plotted and compared once again, on the left those for 3ftAfc, on 
the right those for ^A^. For '^Ai the value of uji is slightly larger for the broader cylinder; 
the other values are very close. For '^Ai the value is smaller for the narrower cylinder, 
and the other values are again very close. In Figure 123 a similar comparison is made of the 



distributions of 3?Ai on the left and of ^Ai on the right. As is to be expected from conditions 



27| , the distribution of '^Ai is shifted to the right. It is more shifted for the narrow cylinder 
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than for the broad. The results encourage the behef in the markovian hypothesis, even 

though it is hard to imagine that experiments as coarse as these could ever successfully 
refute the hypothesis because some shielding is inevitable. The question is rather how much. 
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4 Cylinders of variable length and the phase. 



We have seen in Paragraphs 2]3| and p.2| that the measures m£,,c can be used to recover 
the conformal exponent associated to the spin-spin correlation at the boundary and in the 
interior. Various formulas in the theory of free fields suggest that critical exponents might 
also be obtained from the analogue for the field h of the variable x defined in Paragraph 
for the free boson 0. We refer to this variable as the phase, and our examination in this 
section, although brief, indicates clearly that it also can be used to reproduce exponents of 
the classical Ising model. 

The variable x for the boson field measures the difference between the constant terms 
in 01 and 02, the restrictions of to the two boundaries of the cylinder. It takes its values 
in the interval [0, 2ttR) where R is the radius of compactification. (See Paragraph An 
analogue for the Ising model on the cylindrical LV x LH square lattice Q\j is defined using 

p 

where S is the unit vector in the horizontal direction and the sum runs over all sites p in the 
lattice that have a right neighbor. Because the jumps of h are chosen at random between 
±7r, it is natural to study the distribution of 

X = x' mod 2TT 



instead of x'. The normalization of x is such that a closed curve of discontinuity in h that 
wraps around the cylinder, in other words that is noncontractible, gives a contribution of ivr 
to X. Clusters intersecting the boundary contribute ttA/LV to x where A is the numbers 
of boundary sites inside the cluster. However contractible curves surrounding clusters of 
constant spins not intersecting the boundary do not contribute. 
In Section |^ we introduced, for the cylinder D, the measure 

mDi{ak},{h}) = lim lini m^^({afc}, {6^}) 

N—>oo 

defined on the space S^i with coordinates {{ak},{bk}), A; G Z \ {0}. As we observed in 
Paragraph |3]^, this can also be regarded as a measure m£)(-0i, -02) on a space of distributions, 
one ipi on the circle at one end of the cylinder and one -02 on the circle at the other end. We 
could as well have defined 

mD{i'i,i^2,x) = mD{{ak},{bk},x) = lim Mm m°[f {{ak} , {bk} , x) 

taking the variable x into account. The probability ?72,£i(0'i, -02) is a conditional probability, 
thus - speaking imprecisely - we have integrated over the variable x. Writing all measures 
informally as measures absolutely continuous with respect to a Lebesgue measure on the 
underlying spaces, we express this as 



dmoiiak}, {bk}) = dniDii^i, ^2) = ZdH^i, ^2) #1 #2 
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with 



Zd{^Pi,iP2) = / Zd{iPi,^P2,x) dx. 







This is a convenient notation that avoids technical explanations about conditional probabil- 
ities and also reminds us of the connection between the measures and partition functions. 

4.1 The measure mq{x). 

We first consider dmq{x) = drriD^x) = Zy{x)dx, D = D{q), with 

Zg{x) = J ZDi{ak},{bk},x)Yldakdbk = j ZD{^i,tp2,x) #i d^2, 

the choice between the three notations ({o-jt}, {bk}), ({^fc = ikak}, {Bk = ikbk}) and (V^i, 
being a matter of convenience. We shall parametrize by the variable q the cylinder D in the 
plane of length I A, I = ln(l/g), and circumference 27^ A, with A arbitrary. It is mapped to 
an annulus the ratio of whose inner and outer radii is g by z — > exp{z/A). The measure is 
normalized 



Zq{x) dx 



and its Fourier expansion is 



ikx 

2ti 

k^O 

We can try to expand each coefficient in a series of powers of q 

00 

j=0 

We expect from the original calculations on the Ising model or from arguments of conformal 
field theory that ao = 0, although we admit Ck{ao) = 0, and that ai = |. The remaining 
Uj should be at least |. (The usual argument of conformal field theory would select the 
exponents 0, | and 1, and all those differing from these by two positive units, but it requires 
unitarity. It is not yet clear to us to what extent unitarity is pertinent in the present context. 
The whole Kac spectrum could intervene - at least our experiments are not fine enough to 
rule out a = I which is smaller than a = 1.) 

We have run two sets of experiments to measure the smallest exponent in z/i(g), one for 
LV = 59, the other for LV = 117. As g — > 0, that is for long cylinders, the graph of Zg 
is practically of period vr, instead of 27r, and the odd Fourier coefficients C2fc+i(g) — > 0. 

The physical reason for this behavior is that, for very long cylinders, several noncontractible 
curves of jumps in h are likely to occur and configurations with an even or an odd number 



of these curves will arise in approximately the same numbers. Figure |25| shows, for the long 
cylinder of size 117x801, the distribution of the variable x' (before the identification x' ~ x' + 
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Figure 25: The distributions of the variables x' and x for the cyhnder 117 x 801. 

27r) and of the variable x. The peaks for x' are centered on the integer multiples of vr, clearly 
underlining the role of noncontractible curves of jumps. The figure shows configurations 
with n curves, |?t,| = 0, 1, 2, 3, 4, and the data also indicate that |n| = 5 and 6 were obtained 
in the sample of 1.6 x 10^ configurations. Even for |n| = 4 the probability is fairly large. 
It should be remembered that only ^ of the configurations with 4 noncontractible curves 
will contribute to the peak around 47r. The distribution mq{x) is, for this cylinder, almost 
perfectly periodic of period vr. 



Figure 26 is a log-log plot of ui as a function of q. The data for the cylinders with 
LV = 59 are marked by "•" and those with LV = 117 by "+". The shortest cylinders were 
59 X 27 and 117 x 53. We measured several other longer cylinders for both LV^s. We decided 
to discard for both the figure and the fits the measurements of ui whose 95% confidence 
interval was more than 5% of the measurement itself.0 The linear fits of the log-log pairs 
give a slope of 0.12506 for LV = 59 and of 0.12478 for LV = 117. The line on the figure is 
the latter fit. The value ai = | appears clearly. We did not check its universality but there 
is no reason to doubt it. 

4.2 The ratio Z+-{q)/Z++{q). 

Let Z^^{q) and Z^^{q) be the relative probabilities that with constant boundary conditions 
on a cylinder of parameter q the spins are equal at opposite ends or unequal. There is a 



•^The Fourier coefRcients are given by — X]?=i ^iPi where Ci — cos((i — ^)2TTk/2LV) and pi are the 
frequencies for the 2LV bins in which the data are distributed. We use Ui for the number of data in the 
i-th bin and N for the sample size. Hence pi = Ui/N. Since the distribution of the rii is a multinomial 
MULT(iV; ni, 722, • • • , ?^2Ly-l), the first moments are (rii) = Npi and {niUj) — N{N — l)piPj + NpiSij. 
Therefore Var(i'fe) = jf^J2i^j PiPji^i ~ ^j)^. For the cyhnder 117 x 801 discussed above {q « —43.0), the 
measured i>i with the 95% confidence interval is 0.00468 ± 0.00110 even though the sample was larger than 
1.6 X 10^. It was not used for the fit. 
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Figure 26: Log- log plot of the Fourier coefficient ui as a function of q. 
well-known formula due to Cardy [|Cl|] , 

Z Xi{q) + V2x2{q) 

with 

Xiiq)= n (1 + ^^")' 

m>0, m odd 
m>0, m even 

We could, in experiments, fix the spins along one or both of the two ends of the cylinder 
to be constant. This leads to alternative measures mq{{bk}, x), in which the spins at the left 
end are taken to be -f-l, and (Z+_(g), Z^+{q)). The question arises whether 

mq{{bk}, x) = mg{{ak = 0}, {bk}, x) (29) 

and whether 



Z++{q)5o + Z+^{q)5n = mg{{ak = 0}, {bk = 0},x). (30) 

These two equations require some explanation. The measure mg{{ak = 0}, {6^ = 0},x) 
is understood, in so far as it can be assumed to exist, to be the conditional probability 
defined by the probability measure mg{{ak}, {bk},x), the conditions being ak = bk = 0, or 
equivalently = 5^ = 0, V/c G Z \ {0}. Experimentally this means that it is a distribution 
that we approximate just as we approximate mg{{ak}, {bk},x) itself except that we discard 
all samples for which the restrictions hi and /i2 at the ends of the cylinder do not lie in 
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a suitably chosen neighborhood of 0. The neighborhood is thus to be as small as possible 
but large enough that we do not reject so many samples that the number of useful samples 
becomes impossibly small. We define mq{{ak = 0},{bk},x) in the same manner, but the 
condition is now that = 0, V/c e Z \ {0}. 
If (^0]) is valid the distribution defined by 

f f Z,{{Ak},{Bk},x)lldAkdBk (31) 

with sufficiently small c^'s should be approximately a(g)5o + b{q)5.„, thus a sum of two 5- 
functions with coefficients whose ratio b/a is given by (PSp. Similarly the distribution 



[ Zq{{Ak},x)l[dAk (32) 

J\Ak\<Ck k 

provides another ratio 6/a to be compared with (pSl). 

Measuring these two ratios b/a is difficult. The ratio Z^_/Z^^ decreases from 1 at g = 
to at g = 1. Large ratios Z+_/Z++, those easier to measure, correspond therefore to long 
cylinders. For these the variables Ak and Bk are independent and their distributions are 
known from previous sections. The effect of the constraints can therefore be estimated by 
using Tfc = Probg=o(|v4fc| < c^). Even by imposing restrictions \Ak\ < Ck and \Bk\ < Ck only 
for k = 1,2,3, leaving the other variables free, a choice of ri = r2 = ~ 0.1 cuts the 
number of admissible configurations by a factor of one million for the measurement of (|3l| ) 
and the measurement is impracticable. For shorter cylinders (g 1), the ratio Zj^_/Z^^ 
drops quickly. For a circumference four times the length, the ratio is less than again 
difficult to measure. We limited ourselves to a small window of r = 1/2it, choosing six values 
corresponding to values of q increasing by a factor of approximately 4 at each step. Table 
IV gives the values of r, q, the (rather small) lattices we used and Cardy's prediction. The 
ratios b/a were measured for the constraints: 

ci ~ 0.377 C2 ~ 0.653 C3 ~ 0.929, (33) 

the others being infinite. These numbers correspond to the following probabilities 

Prob(|v4i| < ci) = 0.2 Probd^a] < C2) = 0.3 Probd^s] < C3) = 0.4 

if the cylinder were of size 79 x 157 like the one used in Section ^ For this long cylinder and 
these constraints applied at both extremities, only a fraction (0.2 x 0.3 x 0.4)^ ~ 0.0006 of the 
configurations would be used. We observed that for the shorter cylinders of Table IV more 
configurations passed the test. The difficulty of getting proper samples for the measurement 
of ( ^21) is of course less acute. 

Three sets of measurements were taken. For the first set the constraints given by ( P3| ) 
were applied at both extremities of the cylinders and is thus of the form (|3TD. In Table IV 
it is refered to as const/const for "constrained". For the second they were applied at one 
extremity while the spins at the other were forced to be the same though they were allowed 
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to flip simultaneously during the Swendsen-Wang upgrades. This corresponds to (|32|) and 
is refered to as const/ fixed. The last set is the measurement of the ratio Z^_/Z^^, that is 
the case fixed/fixed. For each lattice enough configurations (> 20 million in each case) were 
generated so that at least 30000 contributed to the integral (|3T|) . Far larger samples were 
obtained for the two other sets. 





LV X LH 


79 X 122 


79 X 104 


79 X 86 


79 X 68 


79 X 52 


79 X 34 




r = LH/LV 


1.544 


1.316 


1.089 


0.861 


0.658 


0.430 




Q 


0.0000611 


0.000256 


0.00107 


0.00448 


0.0160 


0.669 






0.408 


0.331 


0.249 


0.165 


0.0927 


0.0260 




const/const 


0.419 


0.341 


0.276 


0.193 


0.117 


0.0419 


b/a 


const/fixed 


0.411 


0.338 


0.259 


0.179 


0.101 


0.0301 




fixed/fixed 


0.4071 


0.3289 


0.2494 


0.1640 


0.0916 


0.02539 



Table IV. Ratios b/a measured for several cylinders. 

Because of the small sample, especially in the case (|3l|) , large statistical variations are 
expected between neighboring bins and smoothing provides an efficient tool to identify the 
two local maxima around x = n and x = whose ratio was used as a measurement of 
b/a. These measurements appear in the last lines of Table IV. (Smoothing was done as 
in Paragraph |2.2| . The smoothing parameter was chosen as if the distribution of x were 
approximately the sum of two gaussians centered at ^ = and 9 = tt. The ratios b/a did 
not seem to be very sensitive to the exact choice of the smoothing parameter. Of course the 
case fixed/fixed does not require any smoothing since the distribution is actually of the form 
0'{q)^o + b{q)6Tr.) The measurements for constrained/ constrained and constrained/ fixed are 
systematically larger than the predicted values though they are very close, in fact closer for 
longer cylinders than for shorter ones. 

It is useful to see how the choice of constraints changes the measured ratios b/a and 
whether the distribution of the variable x is at all similar to the proposed sum a6o + M^r- 



For the cylinder 79 x 52 we compared four sets of constraints for the measurement of (31). 
The first set consisted of no constraint at all, that is all the c^'s were infinite. The second 
was the one used before and the finite CkS for the third and fourth sets were 



Cl ~ 0.259 C2 ~ 0.653 cg ~ 0.929 



and 



Cl ~ 0.259 C2 ~ 0.441 cg ~ 0.614 C4 ~ 0.782. 



These c^'s correspond to ri = 0.1, r2 = 0.3, ra = 0.4 and ri = 0.1, r2 = 0.15, = 0.2, r4 = 
0.25. For the fourth set only 3152 configurations were admissible out of the 200 millions 
generated and they were distributed in the 2LV = 632 bins. Errors are large in this case. 
Instead of smoothing as before we compared the four sets by expanding their histograms in 
Fourier series keeping only the first ten terms. The ratios b/a are sensitive to the number of 
terms kept. Only the first two digits of the ratios given below, at the end of this paragraph, 
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Figure 27: The distribution of x for four different sets of constraints {ck}- 

are reliable. The smoothed distributions are shown on Figure If the distribution goes 
to aSo + as the constraints become more stringent then the peaks at and vr should be 
narrower and the distribution around | and ^ should go to zero as one goes to the first to 
the fourth set. This is what happens with the four curves. At | the top curve is that with 
no constraint and the one closer to zero corresponds to the fourth set of constraints. Even 
though the values of a and b for the three last sets are quite different, as they should be, 
their ratios are strikingly close: 0.120,0.109 and 0.114. 

Finally we compared the ratios b/a for the three lattices 79 x 52, 158 x 104 and 316 x 208 
using always the constraints (pSl). The numbers of admissible configurations were 40409, 
9931 and 8816 and the ratios b/a, obtained again after truncation of their Fourier series, are 
0.120, 0.129 and 0.130. These numbers are the same within the statistical errors though the 



values of a and b are again different. Figure 28 shows the three distributions, the sharper 



peaks being for the smaller lattices. It seems that smaller c^'s are necessary for finer lattices 
if the peaks are to be as sharp as for the coarse lattice. 

It is not clear whether the above measurement technique can reproduce accurately the 
ratios with a proper choice of the c^'s and the size of the lattice. The very 

superficial analysis we have done does not indicate any decrease in the small gap appearing 
in Table IV for the short cylinders. Still the measurements and the predictions are very 
close. 

4.3 The measure mq{{ak} , {bk} , x) for long cylinders. 

Some identities are suggested by the previous experiments. For infinitely long cylinders the 
following hypothesis seems natural 

Zg=o{{ak},{h},x) = Zg=o{{a_k},y - x)Zg=o{{bk},y)dy. (34) 

Jo 
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Figure 28: The distribution of x for the constraints (^) on the three lattices 79 x 52, 158 x 104 
and 316 x 208. 



As evidence, integrate with respect to the and bk- On the left we obtain 

z/fc exp{ikx) 



and on the right 



if 



271 l/ifcp exp{ikx), 



k 

We have, by definition, z/q = /io = 2^ and z/^ = /i^ = if A; is odd. Experiments on a cylinder 
with 59 X 401 sites yield 

z/2 ~ .00273 fX2 ~ .0208 27ifxl ~ .00271, (35) 

z/4 ~ .0000267 /i4 ~ -00279 2nfxl ~ .0000488. (36) 

Unfortunately only the first line carries any conviction. It may not be possible to measure 
z/4 with any accuracy. 

The measure dmQ{{ak} , x) = ZQ{{ak},x)dx may be of some interest, but we cannot offer 
any precise hypotheses. It can be expanded in a Fourier series. 



Zo{{ak},x) = ^iJ,j{{ak})exp{ij 



x) 



in which /io({afc}) = 1 and fij{{ak}) = for j odd. Then, for example, /i2({afc}) is a function 
of {ttk}, or equivalently, of {^fc}, but, in spite of considerable effort, we have no idea what 
this function might be. 
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Figure 29: The second Fourier coefficient function of \Ai 



A simpler function is 

/ ^^2{{ak})lldak = f{\A,\). (37) 



k>2 



The experiments indicate that 



, sin(6x7r) , , 



with a ~ 0.415 and b ~ 0.603, but this can be no more than an approximation, as Figure ^ 
indicates. (It was obtained for the cyhnder 157 x 1067 with a sample of more than a million 
configurations. The error bars are indicated.) 

The functions iJ,j{{ak}) possess little symmetry. They are invariant under a rotation, 
thus under a simultaneous transformation of all variables — ^ e^^^au, 9 arbitrary, but not 
obviously under anything else, so that for example, 

y"2({afc}) JJc^Ofc 

fc>3 

is a function of three variables, |^i|, ond axg{A\/ A2). The functions /ij({aA;}) are 
intriguing, and we would have very much liked to discover more about them. 
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5 Crossings. 



5.1 Events and the two hypotheses. 

Crossings are one of the main order parameters for percolation models. Consider, for exam- 
ple, a rectangle covered by a regular lattice. A configuration is fixed when each vertex has 
been declared open or closed and this configuration has a crossing if it is possible to move 
on open sites joined by lattice bonds from the left side of the rectangle to the right one. A 
probability is usually defined on the set of all configurations by fixing the probability p that 
a site is open, so that a site is then, of course, closed with probability 1 —p. The probabilities 
for each site are independent but equal. In the limit of mesh length zero, the probability of 
such horizontal crossing is known (rigorously) to have a singular behavior as a function of 
p, being for p < pc and 1 for p > Pc, for a certain constant Pc G (0, 1) that depends on 
the lattice. This definition can be extended readily to the Ising model by replacing crossings 
on open sites by crossings on spins of a given sign, say, for example, of positive sign. The 
probability of crossings on clusters of + spins is not a familiar order parameter for the Ising 
model, and it is not even clear that it is not trivial, thus identically 1 or identically 0. We 
examined it, at first, only out of idle curiosity, following a suggestion of Haru Pinson and 
were somewhat astonished to discover that it is far from trivial. With hindsight, it does 
have some immediately appealing features and has been studied before although not with 
the same goals ||KSC|| . It is related to spontaneous magnetization and to the geometry of 



the main cluster. It even turns out to share striking properties of the percolation crossings: 
universality and conformal invariance ||LPS||. Whether a formula for it analogous to that of 



Cardy [|C2|| for percolation remains an open question. We recall the definitions. 

Let D he a. domain and D' a closed subset of D. Let a = {(ai, 02), • • • , («2n-i, «2n)} and 
P = {{Pi, P2), ■ ■ ■ , (/?2m-i, /?2m)} be sets of n and m pairs of intervals in the boundary of D' 
such that the 2(m + n) intervals are pairwise disjoint. (In fact, the intervals need not be in 
the boundary of D' but these are the only cases we treated.) Let a lattice Q be superimposed 
upon the domain D. Let F be a configuration for the Ising model {Q, J) on D and F' its 
restriction to D'. We shall say that the event E specified by the data {D,D',a,f3) occurs 
for the configuration F 

(i) if for every pair {a2i-i, a2i), i = I, . . . ,n, there is a connected cluster of + spins for F' 

that intersects both a2i-i and 0:21 
(a) and if for no pair {P2j-i, P2j), j = I, ■ ■ ■ ,m, is there a connected cluster of + spins for 

F' that intersects both P2j-i and P2j- 
(For percolation the definition of an event is simpler as the introduction of the larger domain 
D is superfluous, so that the measure on the configurations on D' is independent of the choice 
of D. Thus one takes D = D'.) Let Qa be the lattice Q shrunk by the factor a and let vr^''^'*''^ 
be the probability of the event E = {D, D' ,a, (3) for the Ising model {Qa, J) at its critical 

(G 7^ 

point, then vr^ ' will be defined as 

vrf'^) = lim7rf'^)''^ 

a— +0 

if the limit exists. The two hypotheses of universality and conformal invariance are then 
identical to those proposed in ||LPl5| ] for percolation. 
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Hypothesis of universality: For any pair of Ising models {G,J) and {Q',J'), there 
exists an element g of GL{2,M.) such that 



7rf "'^ = for all events E. (39) 

Hypothesis of conformal invariance: Let (^n, Jo) be the Ising model on the square 
lattice with critical coupling J\j. Let (p be a map satisfying the same requirements as in the 
hypothesis of conformal invariance of Section |^. Then 

T^iSnJa) ^ for all events E. (40) 

It is best to observe explicitly that the map acts on both D and D' , so that if D is 
the whole plane there are very few admissible 0. The following two paragraphs describe 
simulations done to examine these hypotheses when D = D' (Paragraph |5.2[ ) or D' <^ D 
(Paragraph p.3| ). 



5.2 D = D'. 

For the first events to be considered we take D = D'. Their description is simple when the 
geometry of D is that of a rectangle. We introduce the notation 7r/i(r) and vr„(r), instead 
of He, for events E occuring on D, a rectangle with aspect ratio r = width/height, with 
a single pair (ai,a2) and an empty (3. For the probability of horizontal crossings tt/j the 
two intervals ai and a2 are the left and right sides and for the probability of vertical 
crossings, the top and bottom. The probability 7r/i„(r) will give an example of an event with 
two pairs a = {{ai,a2), (q!3,Q!4)}. It is the probability of having simultaneously horizontal 
and vertical crossings in a rectangle D of aspect ratio r. Note that the number 7Th{r) — vr/i„(r) 
is the probability to have a horizontal crossing without having a vertical one. It thus provides 
an example of event E with one pair a and one pair p. Finally we introduce vr^(T') and vr;f (r) 
whose corresponding events have a single pair {ai,a2). For vr^, ai is the vertical segment 
splitting the rectangle in two parts of equal areas and a2 the right side. For vr^, ai is 
the horizontal segment in the middle of the rectangle and a2 the top side. For these two 
probabilities, we could also have taken D' to be the half-rectangle bounded by {ai,a2} 
because a path joining a2 to ai reaches ai before it leaves this half-rectangle, so that the 
sites outside the half-rectangle are superfluous. 

Two difficulties limit the precision of the numerical measurements. The first one is the 



limitation due to a choice of convention and was discussed at length in [|LPPS|] . Since vr^^'"^^ 



are approximated by measurements on finite lattices, the exact position of the domain D 
with respect to the lattice must be specified by convention; or, equivalently, a prescription 
must be given for calculating r for a rectangle with LH sites in the horizontal direction and 
LV in the vertical one. To examine the sensitivity to convention consider an extreme case. 
Suppose that in convention I the width is that of the narrowest rectangle containing the LH 
horizontal sites and that in convention II, the width is that of the widest. For the square 
lattice oriented so that its bonds are parallel to the sides of the rectangles, the difference 
between the two widths is 2 mesh units. If both conventions measure the height in the same 
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way, the discrepancy for tt^j between the two conventions is 




the prime denoting a derivative. These numbers can be estimated from the data of Table VII. 
Table V gives an order of magnitude for this limitation on precision for the two probabilities 
TT/i and TT^ at the center (r = 1) and at the extremities (r = 0.1361 and 7.353) of the range 
of the aspect ratio we measured. Our conventions are given in the appendix; whatever they 
are, the above limitation is unavoidable. 



r 




^Itt' I 

Lv\^h\ 


Statistical error 






Statistical error 


7.3 


0.02 


3 X 10-4 


6 X 10-4 


0.12 


1 X 10-4 


1 X 10-3 


1.0 


0.50 


3 X 10'3 


2 X 10-3 


0.66 


3 X 10-3 


2 X 10-3 


0.14 


0.98 


2 X 10^3 


6 X 10-4 


0.99 


5 X 10-3 


3 X 10-4 



Table V: Sensitivity to convention and statistical errors for a sample of 200000 for 3 values 
of r on a lattice containing fa 40000 sites. 

To confirm the conformal invariance we also measure all these probabilities for comparable 
geometries on the disk and the cylinder. The Schwarz-Christoffel map can be chosen so that 
the four vertices of the rectangle of aspect ratio r correspond to the four points ±e^'^ for 
some 9 e [0, f], on the unit circle. Notice that r = corresponds to ^ = |, r = 1 to 
6 = ^ and r = oc to = 0. The slope of the function 9{r) at r = is zero. This means 
that the sensitivity to convention is magnified for values of 9 close to 7r/2. For example 
we measured the probabilities tt in the rectangular geometry for five different values of r in 
the range [0.1361,0.1647] . The corresponding range of 9 is [1.57051,1.57076] and, on the 
disk of radius r = 300.2 mesh units that we used, at most one site can be contained along 
the boundary in this interval. This is even worse for the corresponding geometry on the 
cylinder of size 397 x 793 where the tt's have also been measured. Such measurements are 
too imprecise to be useful and we measured the probabilities, on the disk and the cylinder, 
only for the ^'s corresponding to the forty-one values in the middle of the eighty-one we used 
for the rectangular geometry. The arc between the two smallest as well as the two largest ^'s 
among these 41 values is about 3.7 mesh units. Since we have taken the sites in the angles 
[71 — 9,71 + 9), {~9, 9) to define the pair of intervals (cti, q;2), it is clear that a rather large 
systematic error is to be expected. 

Finite-size effects are the origin of the second difficulty. Fortunately the relation 7r/i(r) -|- 
7r^(r) = 1 is verified for the triangular lattice, even for finite ones. This is a well-known 
identity for percolation and the argument for its validity here is the same. For the other 
pairs {Q,J), this relation is not verified for finite lattices, that is tt*^^'"'-''" + rci^''^'^'"' ^ 1. 
Nevertheless, if universality holds, it should be satisfied for the other pairs {Q, J) in the limit 
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Figure 30: Log-log plot of 1 — tt/j — tt^ measured on a square as a function of the linear 
number of sites. 



of zero mesh. Departure from zero of the quantity 1 1 — ir^^^'^^^ — ni^^'^^^ \ for r = jy is 
therefore a measure of finite-size effects. Interpreted differently, this quantity is a measure 
of the error made on 7r/i(r) when the number vrj^^^'-^^^ is used in its stead. A verification for 

a square domain covered by the square lattice indicates that 



is likely to hold when the number of sites goes to infinity. The log-log plot of Figure 30 



shows that (1 — tt/, — vr^,) and LH are related by a power law. (The five points correspond to 
squares with 25, 50, 100, 200 and 400 sites along their edges.) The slope is 0.437 and unlikely 
to be universal. The crossing probabilities tt/j and iiy on the square lattice for the square 
(r = 1) of size 200 x 200 were measured to be 0.4963 and 0.4964. The gap is of order of 
3.5 X 10"'^, comparable to the value of 27r^(r)/LK at this point. Note finally that, even 
though 7r/i(r) + 7rt,(r) = 1 holds for finite rectangular subsets of the triangular lattice, it does 
not follow that TTh{r) is equal to irl^^'^^^ as finite-size effects could alter both nl^^'^^^ and 
^(LH,LV) ^Yiile keeping their sum equal to 1. 

The five plots in Figures ^ to ^ show all the data available: the probabilities for 81 
values of the aspect ratio for the rectangles and 41 for the disk and the cylinder. For the 
rectangles, 4 different Ising models were studied: the three regular lattices with isotropic 
coupling and the square lattice with the anisotropic coupling used in Section ^. Each figure 
contains therefore six sets of measurements, four for the rectangles, one for the disk and one 
for the cylinder. The cylinder is treated as though it were infinitely long and the crossings are 
from an interval on one end to another disjoint interval of the same length on the same end, 
the intervals being chosen so that their position on the cylinder is conformally equivalent 
to that of two opposite sides on a rectangle. Because of the large amount of information 
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Figure 31: log7r/i/(l — TVh) and log7r^/(l — n^) as a function of logr. 



on these figures, the error bars were not drawn. Some of the 95% confidence intervals for 
the measurements of Tih were hsted on Table V and the difference of the extreme values of 
these intervals is equal to 0.07 at r = 0.136, 0.02 at r = 1.000 and 0.07 at r = 7.351 for the 
variable log 71^/(1 — Hh) that appears in Figure |31|. (For the square lattice, the confidence 
intervals on the probabilities are a factor smaller since the sample was 5 times larger.) 
The vertical dimension of the dots on this figure is approximately 0.065 and thus comparable 
to the statistical errors or larger than them. 

In all the figures, one sees clearly some spreading of the data at the two extremities of the 
range of r. The data for the disk and for the cylinder also fall slightly beside those for the 
rectangles around the extreme values of their range (logr ~ ±1). These small discrepancies 
can all be explained by the above two limitations. First, for all the pairs {Q, J) but the 
isotropic triangular lattice, the quantity 7Ch + tt^ is less than one. It is thus likely that finite- 
size effects tend to decrease both TCh and tt^. Since for log r ~ ±2, one of the linear dimensions 
of the rectangle is half what it is around logr ~ 0, the values of tt^ and vr^ should be spread 
more at the extremities than at the center of the range of r; and ir^f^^ should be the largest 
of all measurements. This is what is observed though the spread is noticeable only when the 
small linear dimension is in the direction of the crossing. Second, by keeping the sites inside 
the sector {—6{r),6{r)) or (vr — 6{r),7r + 0{r)), the number of sites (necessarily integral) is 
underestimated, leading to probabilities lower than what universality would predict. This is 
again what is observed. But these discrepancies are rather small. As can be seen from the 
figures the agreement is remarkable. 

Only for the isotropic Ising model on the square lattice is 7]-^^^'^^^ strictly equal to 
^{LV,LH) ^ It is then sufficient to measure the five probabilities tt/j, tt^, tt/i^, tt^, tt^ for 41 values 
of r to cover the same range. We profited from this coincidence and substantially increased 
the sample in order to measure the probabilities with very high accuracy. In this case each 
sample contained at least one million configurations. For the other models we used samples 
of at least 200000 configurations. As can be seen from Table V, even the smaller sample 
size yields statistical errors at worst of the same order of magnitude as the sensitivity to 
conventions. Table VII lists the crossing probabilities tt/j, vr^,, 7r/i„, tt^, tt^ for the isotropic 
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Figure 32: log function of logr. 




Figure 33: log7r^/(l — tt^) and log7r^/(l — tt^) as functions of logr. 



59 



Figure 34: Fits of the asymptotic behavior of tt/j: (a) log7r/i(r) and (b) log(l — nh{r)). 



Ising model on the square lattice; Table VIII lists them for the triangular lattice. This table 
gives an idea of both the difference between the various probabilities as measured for two 
different Ising models and the isotropy of the probabilities: the pairs (tt/j, tt^) and (vr^, vr^) 
are approximately symmetric under the exchange of r ^ even though the lattice is not 
invariant under a rotation of ^. 

For percolation, Cardy's formula predicts the following asymptotic behavior 

log 7r)J'^'^™(r) — > — —r + constant, 

r— >oo 3 

or equivalently 

71 

log(l — vrF'^'^'^°(r)) — > h constant. 

r^o 3r 

The data for the Ising model behave similarly. We used those for the triangular lattice since 
they respect closely the relation 7ih{r) + Tc^^r) = 1. We rejected the ten points at both 
extremities of the spectrum of r because they carry the largest finite-size effect. The 30 
remaining points with largest r were fitted to log7r/j(r) ^ a + br and the 30 with smallest r 



were fitted to log(l — 7r/i(r)) c + d/r. The fits appear in Figure ^ The constants h and d 



turned out to be — 0.16727r and — 0.16647r. A natural guess for both constants is — 7r/6. 



5.3 D' ^ D. 

We measured the crossing probabilities from one curve Cj to another one Cj on the cylinder, 
< « < J < 4, and for the corresponding configurations on the disk. (The curves Ci have 
been introduced in Section |^.) The simulations were done on the cylinder with 397 x 793 
sites and on the disk of radius 300.2 mesh units. The results are tabulated in Table VI. In 
each cell the number on top is the probability for the disk, the one on the bottom that for the 
cylinder and, again, the vertical bar "|" is used as in Section Q to give the statistical errors. 
The agreement is convincing even though the probabilities for the disk are systematically 
larger than those for the cylinders. Again the geometries of the disk and the cylinder are 
not quite conformally equivalent. Only if the cylinder is infinitely long can one hope to have 
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perfect agreement. Since the relative gap increases as the two curves Cj and Cj move closer 
to the middle of the cylinder, the shortness of the cylinder is a likely explanation for the 

discrepancy. 

Since the numbers of Table VI are all close to 1., one more example of crossing probability 
was measured. The event E for the cylinder (= D) is given by the following data: the domain 
D' is delimited by the curve C2 and the right-hand side of the cylinder and ai and 0:2 are 

the two intervals on C2 that correspond to the forty-seventh value of the aspect ratio r 
considered in the previous paragraph (r = 1.35). The data for the disk arc the conformal 
images of those of the cylinder. For the disk and the cylinders the numbers tte are 0.412|1 
and 0.4096|6 respectively. 







C2 


C3 


C4 


Co 


0.999981 1 
0.9999791 5 


0.9976|1 
0.99741|6 


0.9634|4 
0.9631 2 


0.8465|9 
0.84564 


Ci 




0.99958|5 
0.99932|3 


0.9727|4 
0.9699|2 


0.8541|8 
0.85104 


C2 






0.9848|3 
0.9827|2 


0.8643|8 
0.86144 


C3 








0.8995|7 
0.89734 



Table VI: Crossing probabilities from one curve Cj to another Cj for the disk and the cyhnder. 

Another interesting choice is D' ^ D = M.'^. That means measuring crossings on domains 
D' in the bulk. We have seen that the uj^^^^ are larger than those at the boundary by 
approximately a factor of 3. The corresponding variances are consequently smaller and 
the number of large clusters intersecting the central meridians of the cylinder is also smaller. 
Arc there enough of them to break crossings? Or is 7rJ^"^''^(r) a trivial function, namely equal 
to I for all r? 

Such a measurement would amount, in an ideal situation, to thermalizing an infinite 
lattice and then measuring crossings on finite D' inside this lattice. Only the usual 
limitations (convention and finite size of D') would then have to be dealt with. To do the 
actual simulations, the first idea is to truncate D to a finite though large lattice and to 
choose D' as the largest domain possible inside a region in which the behavior of the spins 
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is as close as possible to the bulk behavior. With our present computers, a lattice size of 
practical use contains about 10^ sites. If M? is approximated by a square lattice, then it 
would be of size 1000 x 1000. The domains D' used in Paragraph ^]2| contained around 
40000 sites and the domain D' with r = 1 was therefore 200 x 200. If we compare these 
sizes with disks, as we are interested only in orders of magnitude, the boundary of D' would 
correspond to a circle of radius one fifth that of D. The distribution m£, Q£,i on the boundary 
of D' is approximately equivalent to that of a circle at a distance of 100 mesh units from 
the boundary of the cylinder 397 x 793. Figure |1^ (Paragraph |3.2| ) shows that the first four 
Fourier coefficients are still far from their bulk distribution. These coefficients are precisely 
those measuring the large clusters responsible for creating crossings or for breaking them. 

But as we have seen (Paragraph p.2|) , the middle of a long cylinder provides a better 
approximation to bulk behavior. So we confine our experiments to cylinders. If a square D' 
of size 200 x 200 is located in the middle of a cylinder of size about 397 x 793 as before, its 
distance from the boundary is about 300 mesh units and its spins behave essentially as in 



the bulk as can be seen in Figure |I6[ This choice has one possible drawback. It spoils the 
symmetry between horizontal and vertical directions. The mean width of the largest cluster 
is surely not equal to its mean height on a long cylinder. Fortunately a simple quantity, 
7r/i(r) — vr„(^), can be used to quantify this symmetry breaking. 

To enforce the relation 7r/j(r) + 7r^(r) = 1, we took the measurements on triangular lattices 
with 426 X 737 and 852 x 1475 sites, the 737 and 1475 sites being in the longitudinal direction. 
On these lattices, the crossings tt/i, tt^, and vr^ were measured on rectangles with the 81 aspect 
ratios r used before. To keep the rectangles safely in the bulk, we used domains D' with 
approximately 10000 sites. (We used the same domains on both cylinders. See below.) The 
longest rectangle (r k, 7.3) has 40 x 253 sites and its distance from the boundary, for the 
cylinder of size 426 x 737, is similar to that of the square of size 200 x 200 square in a cylinder 
of size 397 x 793 discussed above. The highest rectangle (r ^ 0.13) has 293 x 34 sites and its 
height takes up more than | of the circumference of the smaller lattice, possibly too large 
a fraction if the symmetry breaking is important. The larger lattice helps to address this 
question. We also measure the crossings vr^, and inside a disk of radius 100.2 whose 
center is within one mesh unit from the central meridian of the cylinders. Note that the 
hypothesis of conformal invariance stated above does not relate the crossings in the bulk on 
the rectangles and on the disk. As emphasized, the map must act on both D and D' and 
there is no conformal map from the plane {D) to the plane taking a rectangle {D') to a disk. 

Figures and ^ present the results. Squares (□) were used for the crossings on rect- 
angles and circles (o) for those on the disk. White symbols are for the 426 x 737 lattice and 
black for the 852 x 1475. The two samples were 895000 for the 426 x 737 cylinder and and 
227000 for the 852 x 1475. Even though these data look almost identical to those presented 
in Paragraph |5.2| (Figures 31 and 32), the vertical scale is different. When D = D', tt^ ranges 
from 0.02 to 0.98 as r decreases from 7.3 to 0.14. Here, in the bulk, iih goes from 0.23 to 
0.76 for the same interval of r. 

There is a definite breaking of the horizontal-vertical symmetry. The graph of nhv for 
the cylinder with 426 x 737 sites is clearly asymmetrical. For the rectangles the quantities 
7r/i(r) and 7r^(-) that should be equal if the symmetry was present differ by about 6% for r 
large or small and by 1% for r ^ 1. For the measurements on the disk their departure from 
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Figure 35: logTT/^,/ (1 — tt/^) and logTr^/ (1 — tt^;) as a function of logr. 
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Figure 36: log7r/i„ as a function of logr. 
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Figure 37: The Riemann surfaces Ri and R2 and their corresponding neighborhoods. Only 
one sheet of Ri is presented. 

symmetry varies from 3% to 7%. For both cases, rectangles and disk, the vertical crossings 
are always larger than the corresponding horizontal ones. Large clusters wrapping around 
the circumference are more likely than clusters having about the same number of sites but 
that fail to surround the cylinder simply because the former have fewer peripheral sites than 
the latter. This difference seems to play a role here. If this is so, a better measurement of 
the vr's would therefore be obtained by, say, doubling the linear dimensions of the cylinder 
while keeping the number of sites in the domains D' unchanged. This is why we studied 
the larger 852 x 1475 cylinder. For this new experiment, the asymmetry is essentially gone. 
For example most of the quantities 7f/j(r) and Tivil) differ now by less than 0.5%. Still the 
data for the two lattices remain very close and experiments with smaller cylinders show that 
the curves in Figures ^ and ^ do not change much with lattice size, so that we can assert 
with some confidence that the crossing probabilities in the bulk are well defined as the mesh 
goes to zero, in other words, the crossing probabilities are defined even when D is the whole 
plane. The data, especially those for the 852 x 1475 cylinder, must represent a very good 
approximation to the crossing probabilities in the bulk for the rectangles and the disk. 

In particular, the curves for the crossings on the rectangles and on the disk are now 
distinct and their difference does not seem to be due to the limitation of the experiments. 
There was in fact no reason at all to compare them or to use the parameter r to describe 
the arcs on the disk, for it pertains to a conformal transformation 0^ from the disk to the 
rectangle of aspect ratio r that is no longer pertinent. Nevertheless, it does appear that 
^rectangle < 7rf'^'^(r) for r > 1 and that vr^^ctangiej^^^ > TT^'''^{r) for r < 1, inequalities 
for which we have no very persuasive explanation. For each r, the map (pr extends to a 
conformal equivalence (pr of a double covering Ri{r) of the plane, or rather of the Riemann 
sphere, ramified at four points with a torus R2{r). Figure ^represents Ri and R2 and their 
corresponding neighborhoods. Only one sheet of Ri is depicted here; the other is identical. 
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all data being primed (2 — > 2', etc.). Ri is a double covering of C with cuts tying the four 

singular points on the unit circle uj,ui,—uj and —ui. The cuts were drawn along the unit 
circle. Consequently the neighborhood to the left of the domain 2 on the first sheet is the 
domain 6' on the second. The thicker lines on i?2 are not cuts but circumscribe the images 
of the disks on Ri. The top and bottom sides of the whole rectangle are identified as are the 
left and the right. Both Ri and i?2 are tori. The conformal class of Ri and R2 depends on 
r. The hypothesis of conformal invariance does apply to -j^'^^^^'^^^^^^ g^j^j ^rcctangiccK2(r)^ They 
are expected to be equal. We do not know what relations might subsist between Tj-^ectangie j-^ j 
and 7^^°'=*^°siecii2(0 between ^^^^^'^^^^^^ 7r'^^^^{r). 
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6 Comparison with free fields. 

If we take t] = (pR the interaction for the free field on a square lattice is 

the sum running over all pairs of nearest neighbors. In the continuum limit this becomes 
formally 

/•l/i-\\2 /o\2 



We observe that there is an inconsistency in between the discrete and continuous hamil- 
tonians. For consistency the denominator in (4.3) of that paper has to be replaced by Att. 
We have used the formulas based on the continuous hamiltonian. 

There are at least two properties of free fields that appear again in other models. Either 
might be chosen as a basis of comparison and a means of studying these models. The property 
commonly chosen is the asymptotic behavior of correlation functions. In particular, in the 
plane, 

Mp)-Vm')--^^\PI (42) 
9 

where \p\ is the distance between x and 0; and on a cylinder of circumference 2n 

{ivip)~vm')-ci + -\p\, (43) 

9 

if p and the origin lie on the same generator and |p| is the distance between p and the 
origin with respect to the metric that yields a circumference of 27r. We shall briefly recall 
below the pertinent calculations. In the formula a is a constant that depends on the mesh. 
It could approach infinity as the mesh approaches zero. 

Another property is described in Q. Consider the partition function Z{(f)) with boundary 
conditions, either on a disk so that is a function on the circle, defined however only modulo 
constants, thus for simplicity with constant term 0, or on a cylinder, taken to be infinitely 
long, so that is really a pair of functions 0i, 02, and a constant x, taken modulo 27r. In 
the notation of Paragraph |2TT| 

For the disk, 

Z{<p)=exp{-gJ2\Ck\Vk) (44) 

fc>0 
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Figure 38: The quantity {{h{p) — /i(0))^) — 1.5 In \p\ measured on disks of radii 100 and 300 
covered by a square lattice. 

and for the cyhnder, 

Z(0i,</)2,x) = Z(0i)Z(02). (45) 

Thus, as far as the variable x is concerned, the measure is homogeneous, a behavior that 
constrasts with that of the Ising model discussed in the previous section. 

For the Ising model on a triangular lattice the SOS-model constructed in Section || is 
almost the same as the SOS-model attached, as in for example, to the 0(l)-model on a 
hexagonal lattice,^ and for this model there are familiar arguments that suggest the behavior 
(^) with g = gi = 4/3. We have not tested carefully the universality of the behavior or of 
the constants. Crude experiments for the square and the triangular lattice suggest that the 
behavior is universal but we are not certain that the constants do not vary slightly. 

The function rj of the free-field theory plays the same role as the function h of our 
construction so that to test (Sa) we examine {{h{p) - h{0)y). For what they are worth, the 
results for the plane appear in Figures |3^ and ^ in which the value of 

{{h{p) -h{0)Y) -1.5 In \p\ (46) 

is plotted against \p\. For the square lattice in Figures ^ the experiments are performed 
in disks of radii 100 and 300, an edge of the lattice being taken in each case as unit. The 
experiments are perhaps not to be taken too seriously because the finite size leads to an 
ambiguity. Not only are the states in a disk qualitatively different at the boundary from 
those in the true bulk limit but also the jump lines that in a disk terminate at the boundary 

*In Q the partition function for the 0(l)-model is expressed as a sum over weighted closed curves in 
the hexagonal lattice which is dual to the triangular lattice. Every state of the Ising model leads also to a 
collection of closed curves, formed from the dual edges separating sites of different spin. The weight of the 
collection as a whole can be taken as the mass of the set of Ising states that lead to it. Our prescription leads, 
however, for the individual curves in the collection to different Boltzmann weights than the usual complex 
weights determined locally as in g. For the reasons explained in the following section this does not affect 
the relation (^6|). 
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Figure 39: The quantity {{h{p) — h{0))'^) — 1.5 In \p\ measured on a disk of radius 90 covered 
by a triangular lattice. 



could, in some sense, in the bulk turn and pass once again through the disk, so that working 
in the disk increases the statistical independence. The graphs, in which vertical distances 
are drawn at a much larger scale, suggest that the function is approximately constant except 
close to the origin and near the edge of the disk, where the effect of the boundary manifests 
itself. The constant to which one might imagine the difference (^6]) tending has not yet 
stabilized in the diagrams. There is a difference of about .4 in the minimum of the two 
curves. For comparison, a similar curve for the triangular lattice, obtained once again in a 
small disk of radius 90, an edge of the lattice again being taken as unit, is shown in Figure 
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For the cylinder and the same two lattices, square and triangular, the graphs of {{h{p) — 
h{0)y) appear in Figures HO and ETl We have not used the parameter \p\ in the figures but 
rather the parameter k because it is then easier to explain which part of the curve we used to 
calculate the slope (for the square lattice \p\ = ink/LV, for the triangular \p\ = 87ik/\/3LV). 
None the less the data have been so normalized that if the behavior is, as in equation (^31), 
asymptotically a + b\p\ then the slope of the curves in the figures on their middle, linear 
parts and as functions of k is also b. For the square lattice the cylinder is of circumference 
120 and length 2401 in lattice units; for the triangular of circumference 160 and length 1601. 
The slope for the square lattice is about .460; for the triangular it is about .452, which is 
not a number that we can deduce easily from 4/3. These numbers are close; so universality 
of the slope is strongly suggested. 

To obtain Figures ^ and ^ we construct h as in Sections |] and ^ and use the difference 
between the values of h at points on generators of the cylinder symmetrically placed with 
respect to the central meridian and at a distance of k sites from it. Thus k is necessarily less 
than one-half the length of the cylinder. Since we use all generators there is considerable 
statistical dependence. None the less this yields for a cylinder of size 120 x 2401 a very regular 
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Figure 40: The correlation function {{h{p) — h{0)y) on a cylinder for the square lattice 
120 X 2401. 




graph and if we use that part of it between 100 and 1100 we obtain a fit 37.6275 + .4597fc 
from which the statistically generated values differ by no more than two units at any point 
on this interval, so that the slope should be correct to about two parts in a thousand. (The 
curve is in fact slightly concave and the departure from linearity regular. With a quadratic fit 
and a slightly shorter interval we would do much better with the fit but not with the slope.) 
The experiment repeated on the interval [200, 1000] leads to a slope of .4599 but the same 
conclusions. A similar experiment for a cylinder of size 120 x 1201 yields to a slightly better 
fit and similar conclusions with a slope of .4593. An anisotropic lattice of size 78 x 2401 is 
roughly conformally equivalent to a square lattice of size 120 x 2401. Using the points on the 
interval [200, 1200] we obtain a fit of 39.8942 + .4519/c from which the statistically generated 
values differ by no more than 1.5 units. The difference is again not random but not convex. 
We can again conclude that the slope is correct to about two parts in a thousand. The 
difference between the slopes in the symmetric and the anisotropic cases is .0078. In other 
words, it appears that we obtain the same constant. A triangular lattice of length 1601 
and circumference 160 is conformally equivalent to a square lattice of circumference 120 and 
length a little shorter than 1401. Thus it is long enough. On using that part of Figure ET 



in the interval [100, 700] we obtain a fit 20.6918 + .4515fc that is as good as those for the 
square and anisotropic lattice and suggests, for the same reasons, an error of two parts in a 
thousand. 

We observe finally that the experiments described in Section H, in which the analogue of 
Z(0) is studied, yield the behavior (||) and (|5|) with gs = 2R% = 1.4710 ^ gi. 



For the convenience of the reader, we recall briefly the calculations that lead to (|4^ ) 
and (1^). The average {{i]{p) — 1]{0)Y) is taken with respect to the measure defined by the 
weights 

^^p(^-£j^(viP + S)-v{p)?^ =exp{-{Qr],r])), 5 e {(±1, 0), (0, ±1)}. (47) 



p 

The operator Q is obtained from the relation 



^iviP + S) -V{p)f = XI I -2X?7(p + 5) J r]{p). 

p,5 p \ S / 

We are calculating the second derivative of 

(exp(iaA(?7)) exp {-{Qr],!])) drj/ / exp {-{Qr],!])) drj = -exp [-a'^{Q'^\, A)/4) 



with respect to the parameter a, where A is the linear form rj ^ 7]{p) — //(O) or the function 
Sp — Sq. The second derivative is 

(g-u,A)/2. 

This expression is easier to treat when we pass to Fourier transforms. The two delta 
functions of A can be written as integrals of eigenfunctions of Q. Since the operator Q acts 
on T] so that {Qrj){j)) is g/Sii times 



8r/(p) - 2r/(p + (1, 0)) - 2'n{p - (1, 0)) - 2r^{p + (0, 1)) - 2r^{p - (0, 1)), 
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its eigenfunction e27n(pix+p2j/) corresponds to the eigenvalue g{siv? -kx + sm^ ixy) / -k . Therefore 
\{Q~^\,\) becomes at p = (^1,^2) 



TT 

2^ 



1/2 /•I/2 \ 2ni(pix+p2y) 



1/2J-1/2 sin^ TTx + sin^ Try 



■dxdy, 



or 



27r /"^/^ /"^/^ sin2(7r(pia; +P2I/)) 



1/2 J- 1/2 sin^ vra; + sin^ iry 



dxdy. 



(48) 



The integral outside a circle of small positive radius e about remains bounded as |p| — > 00 
and inside this circle the denominator can be replaced by 7r^(x^ + y^). The result is 



2 r dr ['^^ 

— / — / sin^ (r-u cos 6*) (i6', u = \p\. 
Jo r Jo 



(49) 



The integral of 



is the sum of 



/ — / sin^ru COS e) de = 0{u^ rdr) = 0{l) 
Jo ^ Jo Jo 



and 



dr 



l/u 



2-IT 



sin iru cos 9) d9. 



(50) 



Since sin^ ip = \ — \ cos 2^9, 



2tt 



COSU COS 



and Jq{z) = 0{\z\ ^/^), (|49|) can be replaced by 



TT 



TT 



2nJo{z), 



Inu. 



l/u 



Multiplying by 2/g7i' we obtain (^21). 

For a cylinder we treat a lattice that is periodic in the vertical direction (the pi direction) 
with period A, which for simplicity we take to be even. If p = (0, An), n > 0, the analogue 
of (|Sp is 

27r 1 sin^( 



^ x=":^/2 ^ "^1/2 sin^(7ra;M) + sin^(7r?/) 



Once again we drop terms that remain bounded as A approaches infinity. This yields 



g-K \ A 



1/2 



sin^i^TiAny) 



1/2 



^/2-l ^1/2 



x=l 



siii^^TT Any) 



1/2 xVA2 + y^ 



dy 



(51) 
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We examine the second term of (|5T|) using the identity 

sin^ {it Any) ~ 2 ~ 2 "^^^l^^^^?/)- 

The expression obtained from the term 1/2 on the right is independent of n and on close 
examination is seen to behave hke In A, but that is not pertinent here. Since 

/-Va cos{2n Any) y--^ f^^^ cos{27my) 



which upon integration by parts becomes 

1 sin(27m?/) A/2 ^ 1 j^/'^ 2ysm{2Tmy) ^ 

the second term behaves - independently of A - as 0{l/n). This leaves the first term of 
(^) which is n times 

—^^dy 

1/2 y 

For large An this expression is approximately 27c/g. If, however, we measure the distance 
between p and not in terms of the circumference but in terms of the radius of the cylinder, 
the constant 27i/g is replaced by l/^f as in (^). 

Although we have inferred the relation (|46|) from the corresponding relation for the SOS- 
model associated to the 0(l)-model by the construction of our construction of the 
measure on the set of functions h is much more naive and involves no complex weights. As 
a consequence the measure is no longer gaussian. The relation (^2|), with g = 4/3, applied 
to h suggests that, if it were, the appropriate gaussian would be 

9 rfdfV , f df^^ 




A7T\dx) ^ \dy) ^' 

thus that of (PI). The usual formulas for the expectation of the exponential e^^^^^ of the 
linear function X{h) = h{p) — h{0) then suggest, after renormalization, that the correlation 
function of the spins, thus the expectation of q'^Kv)-iK'^) jg 

g-3/41n{p) _ l/p3/4_ 

The exponent is of course not correct. The explanation is presumably similar to that of 
Section |^. It may be possible, although we have made no attempt to do so, to use the 
functions h to construct in the limit a measure on distributions in the plane and this measure 
may very well share some basic properties with the usual gaussian measure, but it will not 
be gaussian. 

It should perhaps be observed that the random variable — 0(0) is not well defined 
on distributions, so that the expectation of e*'^^^)"*'^*^'') makes no sense. Strictly speaking, one 
should take a smooth function A*^ = q approximating as cr — ^> the difference 5p — 5q of 
two (5-functions, calculate the expectation of e'^(p) of e'^'^^^"^ = e^^''^^\ normalize by dividing 
by the value e'^(po) at a fixed Pq, usually taken at a distance 1 from the origin, and then 
pass to the limit ex — >■ 0. This method was used in Paragraph 
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7 Alternate constructions. 



In this section, we examine briefly other conventions and constructions that we could have 
chosen in our experiments. 

7.1 SOS-model jumps. 

If, as indicated in the introduction, the aim is simply to develop the circle onto the line, 
thereby turning the Ising model into an SOS-model, the particular construction chosen is 
somewhat arbitrary. We could, apparently with equally good reason, replace the jumps of 
±TT by jumps from a set, { — {2k + l)n, —{2k — l)7r, . . . , {2k + l)7r}. A; G N, each choice being 
assigned a probability on which the only conditions are that the sum of the probabilities is 
one and that the probabilities of jumps by equal amounts in opposite directions are equal. 
It is not, at first, clear what effect this has. 

Thomas Spencer pointed out to one of us that the behavior, for jumps of ivr, 

((%)-M0)r)~3/21nb| 

is a consequence of a more geometric hypothesis.0 To construct the function h attached to a 
particular state of the Ising model, we construct curves separating the regions in which the 
spins take different values. Let, in the plane, = N{p) be the number of curves separating 
p from the origin. The hypothesis is that 

(iV(p)) -c^lnH. (52) 

Since h{p) is then obtained by assigning independent values to the jumps of ivr, it is clear 
that Cat must be 2/gTi'^. For a cylinder the analogue of (|5^) is 



{N{p))^c\p\. (53) 
Once again, out of curiosity, we tested this hypothesis numerically for the square lattice. 



The results are presented in Figure ^ in which {N{p))/ In \p\ is plotted for the square lattice 
and two disks of radii 200 and 300. It appears that except at the center and near the 
boundary the quotient is approximately constant but that it is only very approximately 
equal to 3/27r^ ~ .15199. There are several possible causes - in addition to a departure from 
gaussian behavior. As we saw in Paragraph the bulk state is approached only slowly in 
a disk. Moreover the finite-size effects that appear in the examination of (^) appear here 
too. The first consequence is that there will be a tendency to overestimate the number N{p) 
when \p\ is not small in comparison with the radius because the curves in a disk that reach 
the boundary are not allowed to close. In principle, this effect should, for a given |j9|, be 
mitigated as the radius grows. On the other hand, rather than increasing toward .15 as we 



'''The hypothesis (|5^) refers only to the weights attached to curves without regard to orientation and for 
them our weights are the usual ones. When deducing (^6|) from ( |5^ ) the relative weights, complex or not, 
attached to the two possible orientations are irrelevant. All that matters is that they be independent from 
curve to curve. 



75 



0.2 
0.18 
0.16 
.14 
0.12 



0.2 

0.18 
0.16 
.14 
0.12 



50 100 150 200 



50 100150200250300 



Figure 42: The quantity (A^(p))/ln \p\ for the square lattice on disks of radii 200 and 300. 



pass from a radius of 200 to one of 300, the minimum of the curve, decreases from about 
.14 to about .13. Since the smallest pertinent value of \p\ is about 75 and ln(75) ~ 4.3 
and the difference in (]46| ) does not, as we saw in the previous chapter, approach a limiting 
value rapidly, if it approaches one at all, a decrease in the minimum of .4/4.37r^ ~ .01 is not 
completely unreasonable. No conclusions are possible without further study. Our purpose 
here is not, however, to examine (0) but rather to acquire a rough understanding of what 
we might have discovered if we had chosen the jumps in a different way. 

The advantage of ( |5^ and (^3p is that they make clear that the behavior (^2p and 
does not change when the definition of h is modified. If there are jumps of {2k + l)7r with 
probability Wk, k ^ 1^, then persists with a new constant 



cn 



oo 



Wk{2k + lf 



(54) 



There is a similar change in (^). 

The behavior of the functions Z{ip) appears to be quite different. We have performed 
a few rough experiments, replacing the jumps of ±7r by jumps of —Sir, —tc, tt, 37r, each with 
probability 1/4 and by jumps of — 57r, —Sn, —tt, it, 3tt, 5n, each with probability 1/6. If the 
measures continue to exist, but with qb modified as suggested by ( ^If ) then the Fourier 
coefficients would continue to be distributed as gaussians but with g = qb oi (^) multiplied 
by I and by ^ respectively, so that the ideal value of a/ 2R\ / tt ~ .68 of the distribution 
at would be multiplied by a/1/5 or a/3/35 yielding ~ .31 and ~ .20. In the first row of 
Figure]^ (four and six jumps) the distribution of the Fourier coefficients y{Ai for a cylinder 
of size 299 x 599 is compared in each of these cases with a gaussian with the same value 
at 0. There is some similarity but considerable difference. Moreover the value at is close 
to but different from the suggested value. For the higher coefficients the distribution looks 
more and more like a gaussian. In the second row of Figure |^ the distributions of D^As, 
normalized so that the factor \/k with which we are familiar from Section ^ are compared 
with gaussians. Not only are they closer to gaussians, but the values at are closer to those 
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Figure 43: The distribution as a function of (first line) and ^A^, (second line) with four 
jumps (first column) and six jumps (second column). 



predicted by (|5^). 

On the other hand, the first row of Figure ^ in which the distributions of SHAi for the 
three sizes 99 x 199, 199 x 399, and 299 x 599 are compared in each of the two cases, suggests 



that the limiting measures may none the less exist. So does the second row of Figure ^ 
for ^RA^. We have, however, as yet made no serious effort to decide whether this is so, nor 
whether these measures could be conformally invariant and universal. 

Another possible "natural" choice for the relative weights of the jumps nn, n odd is given 
by the Dirichlet form. Its discretized form ( ^7]) used in Section 6 suggests that the weights 
-072^+1 and wi of having a jumps ±(2/c + l)7r or ivr satisfy W2k+i = wf'^^ . If wi is fixed by 
requiring that ^jg^ ^^2j+i = 1, then Wi = a/2 — 1. Then the constant gs would be multiplied 
by the inverse of X^iez ^2i+i(2« + 1)^ = 3, that is |. 

In addition to ( ^2]) we also examined, following a suggestion of Thomas Spencer, the 
behavior of 

{{N{p)-{N{p))f)/\n\pl 

whose behavior is pertinent when attempting to establish (|i6| ) rigorously, in disks of radii 
200 and 300. Although the results are not relevant to this paper, they are presented, for the 



curious reader, in Figure Once again, the curves are extremely fiat, but there is a drop 



for the larger radius that has to be explained. 
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Figure 44: The distribution as a function of ^Ai (first line) and (second line) with four 
jumps (first column) and six jumps (second column) on the cyhnders 99 x 199, 199 x 399, 
and 299 x 599. 
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Figure 45: The quantity {{N{p) - {N{p))f)/ln \p\ on disks of radii 200 and 300. 
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7.2 The Fortuin-Kasteleyn construction. 



The Fortuin-Kasteleyn formulation of the Ising model can be used to map the partition 
function of the high-temperature phase of the model to a percolation-like sum over bond 
configurations. To construct the F-K version of an Ising model on a planar graph Q with 
vertices s E S and bonds b E B we shall form the first barycentric subdivision Q' of Q. Thus 
associated to Q are the vertices s, the bonds b, each joining two sites, and the faces /, each 
face / being bounded by sites and vertices. The sites S' of Q' are the sites in S and points 
obtained by choosing arbitrarily from each bond b and each face / a point in its interior. 
Thus, set-theoretically, S' = S U B U F. The bonds B' are pairs consisting of a bond in B 
and one of its ends or a face in F and a bond or vertex on its boundary. In fact, the bonds 
in B' joining a face to a vertex in its boundary are for our purposes superfluous and are not 
included in our constructions. 

The partition function of the original model is taken in the form 

Z = ^Ylexp{J6^^s),a{t))- (55) 

cr B 

Here s and t are the two sites joined by b. Thus, for a square lattice, sinh(J) = 1, J = .881374. 
For a given configuration a, the clusters, in the sense of this paper, are maximal connected 
subsets of S on which a is of constant sign. To obtain a Fortuin-Kasteleyn cluster we remove 
the bonds of these connected clusters one by one with probability l/u, u = exp J (for a square 
lattice 1/z/ = .414214). This replaces the sum (^) by a sum over decompositions of Q into 
subgraphs, each component being provided with a sign. A decomposition is the subgraph 
obtained by keeping all vertices and removing some bonds. 

cr B 

= 5^(1 + (z.-l))'-{(z.- 1)7(1 + (z.-l))^} 

The sum in the second line runs over all decompositions into subgraphs, each component 
being signed, so that a constant spin is assigned to each of its vertices, and so do the sums 
in the remaining lines. From a signed decomposition we can of course reconstruct, from the 
signs alone, the original state of the Ising model. This state has r bonds that join sites with 
the same spin, so that its probability is v"^ . The number of bonds in the subgraph is q and 
the factor (1 — -YQ-Y~^ is the probability that we arrive at it on removing bonds. If we 
now ignore the spins, the final sum in (|5^) becomes 

if c is the number of connected graphs in the decomposition. 
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To construct the function h we associate to a decomposition a state a' on S'. The value of 
a' is 1 at the vertices of ^, at the bonds of Q that belong to the subgraph, but is —1 at all other 
vertices of Q'. Now h can be constructed as before, except that the jumps are to be ±7r/2 
and not ±7r. It turns out to be instructive, at least for the crossing probabilities, to replace 
the probability 1/z/ by a variable probability 1 — yU between and 1. Thus fipK = .585786. 

We have considered only graphs formed by square lattices on either a cylinder (for dis- 
tributions and correlations) or a rectangle (for crossing probabilities). Our aim was not 
to establish conformal invariance and universality for the F-K construction, but rather to 
acquire a provisional understanding of the way the various objects introduced in this paper 
behave under an alternative description of the model. 

There are two ways to define crossings in the F-K construction. If cluster signs are taken 
into account, a crossing is a cluster of sites with positive spins that joins one side of the 
rectangle to the opposite one. The crossing probabilities considered earlier are recovered if 
/i = 1 but the crossing probabilities are zero if = 0. 

If clusters are unsigned, crossings are defined as in bond percolation. This is more in 
the spirit of the F-K formalism and we shall use this definition. Note that both conventions 
are linked: if 7r_|_ and 7r_ are the crossing probabilities over a positively or negatively signed 
cluster, and 7r_(__ the probability that there are spanning clusters of both positive and negative 
sign, then the probability vr that an unsigned cluster crosses is given by the following obvious 
relation, 

TT = 7r_|- + 7r_ — TT^ 

Our crossing probabilities now depend on two variables: the aspect ratio of the rectangle 
r and the probability fi of not removing a link. We studied each of these variables separately, 
varying one and keeping the other fixed. We first took 1 — /i = 1 — fipx = and studied 
the dependence on the aspect ratio. Results for 7ih{r, fipx), the probability of a horizontal 
crossings in Q' on either -|- or — clusters, are shown in Figure The numbers of sites in Q 
inside the rectangles were around 40000 and the samples 250000. The absence of symmetry 
imphes that duality fails. 



The asymptotic behaviour of \og7rh{r, ij,fk), shown on Figure is found to be 



logTChir, fipx) — ^ —0.5027rr + constant, 

a number reasonably close to 7r/2, despite the rather low statistics. 

In the second experiment we measured the dependence on /i of 7ih{l,^). The results 
presented in Figure ^ raise the question whether yU = fipK might be a critical value of the 
function 7rh{r, ^fk) or at least of 7r/i(l, fipx), that is vr/i(l, /x) would be zero for < /i < fipK 
and nonzero for fi > fipK- It is not obvious from the numbers obtained what the limit of 
the function Trh{l, fiFx), fJ' > fJ'FK, is when the mesh goes to zero. 

We examined the correlation functions 



((%)-%))') 
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Figure 46: \ognh{nFK)/ — T^hif^FK)) as a function of logr. 




Figure 47: Fits of the asymptotic behavior of log nh{r, iifk) as a function of r. 
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Figure 48: 7rft(l,/i) as a function of /i for 100 x 100 and 200 x 200 square lattices Q. (The 
curve of the larger lattice is the top one for large 
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Figure 50: The numbers ujk for 1 < /c < 30 for the cyhnders 149 x 339 and 299 x 679. 

both on a cyhnder and on the plane. According to Q one should expect (provided that an 
analogue of (^) is valid) that for correlations in the plane 

((%)-MO)rW~^lnb|. 



This is confirmed by the graphs of Figure 49. In the first the radius of the disk considered 



is relatively small, about 60 bond units; for the second it is 180 bond units. The presence of 
three distinct curves, corresponding to the cases that p is a new site on an old site, an old 
bond, or an old face, while is taken to be a new site on an old, is curious. It appears that 
they remain distinct in the limit of an infinite radius, but their separation remains bounded. 
Once again the scale in the vertical direction is very large; the curves of these diagrams are, 
in fact, extremely fiat except near the ends. 

On the cylinder, the correlation functions behave as 

{{h{p) - h{q)f)^ ~ a(/i) + - g|, 

at least if p and q lie on a common generator. As observed in Section |^, the quantity a(/i) 
is a constant that depends on the mesh and on the nature of the pair {p, g}, on whether p 
or g is a site, bond or face of the graph Q. If the conventions of the equation ( ^31) are used, 
the value of b^fiFx), estimated on a cylinder of size 99 x 699, is close to .26. 

We studied the distribution of the function h for two sizes of cylinder, 149 x 339 and 
299 X 679, but only on the boundaries, not on inner circles. These cylinders are a little short, 
so that about 8/1000 of the samples are such that the sum of the jumps on a circumference 
are not 0, but passing to longer cylinders of size 149 x 449 and 299 x 899, although it reduces 
this fraction to 2/1000 does not change the conclusions. The measures on the boundary 
appear to be gaussian once again, but with a new value of qb that is a little greater than 
3. (Given the behavior of Uk of Figure ^ it is not so clear what is to be. Further study 
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Figure 51: The distribution constructed according to the FK-definition as a function of 
or ^A^. 

might suggest defining it by the asymptotic behavior of Uk.) We plot the values of djk for 
1 < A; < 30 on Figure ^ The results, coarse as they are, are similar to those described in 
Section 0, although there are curious features that advise against hasty conclusions. The 
collection of values for the two cylinders cross at A; = 7. Graphs of the distributions of 'iKAi 
and appear in Figure |5l]. On the left the results for the cylinders of different sizes 
are compared with each other; on the right the results for the largest of the two cylinders 



are compared with gaussians. Figures ^ and ^ together suggest that the behavior of the 



function h constructed according to the FK-definition might have similarities with that of 
the function constructed by the methods of this paper, but we have not examined the matter 
carefully. In particular, we tested neither conformal invariance nor universality. 



7.3 Infinite temperature. 

For the Ising model at infinite temperature, thus for site percolation in which each site is open 
with probability .5, the crossing probabilities cease to be of interest. They are all or, in 
exceptional and trivial cases, 1. On the other hand, the partition functions Z°°{tpi,ilj2,x,q) 



seem to behave much like those at the critical temperature. In Figure ^ we present results 
for the square lattice on cylinders of size 99 x 399 and size 299 x 1199. On the top, the 
results for 9KAi and 91Aio for these cylinders are compared with each other. On the bottom 
the results for the largest of the two cylinders are compared with a gaussian. As in Figure 
|I|, there has been no renormalization of these distributions, so that if the distributions were 
similar to those of Section ^ the ratio of the heights of the two curves would be 1/^0 ^ .32. 
It is about .34, but the cylinders are still fairly small. Although this has no perceptible 
consequences, these cylinders are short enough that about 15/10000 of the sample states 
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Figure 52: The distributions and ^Aiq at infinite temperature. 

yield jumps whose sum along a circumference is not 0, so that the states at the two ends are 
certainly not independent. In Figure ^ the results for the smallest cylinder are compared 
with those for a cylinder on a triangular lattice of size 116 x 401. This is a very stubby 
cylinder, but, curiously enough, once again only about 15/100000 of the states are such that 
the sum of the jumps along a circumference is not 0. Figure ^ is analogous to Figure |^: the 
two sequences of points on the left are for the square lattice, the upper for the smaller of 
the two cylinders, the lower for the larger; the two sequences of points on the right are for 
the smaller of the cylinders with a square lattice (lower set) and for the cylinder with the 
triangular lattice (upper set). If Figure ^ is to be believed the constant 2Rb changes and 
becomes approximately one-half its previous value, but, as with the other examples of this 
section, our aim was more qualitative than quantitative. 

We can also introduce, as in Section ^, the partition functions Z°^_^_ and or the 
measures Z°°{ip,x). A little refiection shows that the two numbers do not depend on q and 
are now both equal to |. The analogue 

^++5o + Z^_6^ = m^({6fc = 0}, x) 

of an equation deduced from equations ( PUj) and (^) is not valid, rather the simulations 
suggest that, if Z°°(0,x, g) exists, then it peaks at tt/2 and 37r/2. On the other hand, 
Z°°(0,0,x, g) has quite a different behavior and a relation between Z'^Q{{ak} , {hk} , x) and 
Z'^Q{{ak},x) like that of equation (|34D is difficult to ascertain because the functions h do 
not very often have level lines that encircle the cylinder, even when the cylinder is very long. 

85 



-7.5-5-2.5 2.5 5 7.5 



-4 -2 



2 4 



Figure 53: The distributions and 3?Aio at infinite temperature compared for square and 
triangular lattices. 
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Figure 54: The numbers oOk at infinite temperature. 
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Figure 55: An unsatisfactory but curious approximation. 
More pertinent to the study of the measures at criticahty is that the behavior of 

■?/'2, x)(i?/'2/ / Z'^{ipi,ip2, x)d'tp2dx = /i^(^) exp(^fca:) 



is similar to that of Zo{iIj, x). This may be of some advantage for numerical studies since at 
infinite temperature no thermalizations are necessary. Consider for example the analogue 



f°° of the function defined by equation (|37|). There is, once again, a simple, rough, but 
inexact - as is clear from Figure - approximation to this function. 



r{x)r^a 



Sm OTTX 



bTCX ' 

but, as before, we were unable to improve upon it in a useful fashion. 
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8 Afterword. 



Most of the phenomena we have examined in the paper are manifestations of the influence 
of the boundary, so that to some extent the thrust of the paper runs counter to the notion 
that statistical mechanics, especially as it refers to critical phenomena, is the study of bulk 
properties. Since critical behavior appears when the appropriate equilibrium between the 
strength of the interactions and the number of paths by which the interaction is transmitted 
over long distances is achieved it is not, from a mathematical standpoint, such a bad idea to 
study criticality by examining the consequences of systematically blocking large numbers of 
these paths. On the other hand, renormalization is usually conceived in bulk terms. So it is 
a relief that the distributions investigated here, whose ultimate purpose is the introduction 
of a concrete notion of fixed point, do not become trivial when the boundary moves off to 
infinity. 

For the Ising model, however, in contrast to percolation or to the free boson, there are 
formal difficulties in the introduction of a closed renormalization transformation that we are 
still unable to overcome, even with measures that continue to have a meaning in the bulk. 

One connection that we would like to make, and that is suggested by some of the ex- 
perimental results of the paper, is with the notion of conformally invariant field theory in 
the strongly geometric form envisaged by Graeme Segal (JS^, |S2l). It may be that the basic 
objects of that theory are constructible from the measures examined in this paper. Recall 
that in that theory one of the first objects to construct is a Hilbert space 

associated to a circle with parametrization. In addition, suppose that we are given a Rie- 
mann surface S with boundary C consisting of disjoint parametrized circles Ci, . . . , Cm and 
C(,...,C^, the parametrizations being given by real analytic functions. Then (equation 
(1.4) of the theory is provided with an operator 

t/s : rr ^ rr. 

We might suppose that Ti. is the L^-space of a measure /i on the space of distributions on 
the parametrized circle. One such measure whose existence is suggested by the experiments 
of this paper is the measure on distributions on a circle in the plane or on the central circle 
of an infinite cylinder described at the end of Paragraph |3.2| , thus the measure defined in 
the bulk. It is possible that L'^^ii) is, if not 7i, then the vacuum sector Ti^y ® Haa or some 
other subspace of TC. 

Consider the annulus of inner radius q and outer radius 1 and the operator Uq associ- 
ated to this surface. Take C2 to be the outer circumference with the natural parametrization, 
Ci the inner, and C to be their union. We consider the annulus as imbedded in the plane 
or, if we treat it as a cylinder of finite length, as being imbedded in a cylinder extending to 
infinity in both directions. The construction of bulk measures suggested in Section ^ yields 
experimentally a measure niq = m^^^c on the product of the spaces of distributions on Ci 
and C2. If, as we might suppose, niq is absolutely continuous with respect to /i x /i then it 

88 



is given by a kernel 



It is not impossible that the operator Uq, or rather its restriction to the sector represented 
by L'^^jj), is given by 



An essential feature of these operators would be the relation Uq = UqJJq,^ when q = qiq2 
which would follow from a relation 



Let C be the circle that separates the annulus of parameter q into annuli of parameters qi 
and q2. We apply the notions of conditional probability and the markovian property, as well 
as the obvious symmetry of Kq, to the bulk measures. Thus 



from which the equation (^) would follow. 

These are tentative suggestions, and we only make them to confess that we have not yet 
had an opportunity to test them experimentally. That may not be an easy matter. Nor 
do we know whether they are confirmed by the conventional wisdom. To construct some 
analogue of Segal's operators on the whole Ti it may be necessary to utilise the phase of 
Section ^, but here again more reflection is necessary. 

Another set of experiments waiting to be performed, although here the outcome is more 
certain and the experiments therefore less tempting, is an examination of the behavior of 
the measures in a neighborhood of the critical point as we vary J (or the temperature) 
and introduce a small magnetic field. The limits as the mesh goes to are expected to 
exist no longer, but the behavior of the measures, of their moments for example, should 
yield the usual critical exponents v and A and should correspond to the usual intuition. 
We are nevertheless curious to see how the geometry of the fixed point is reflected in the 
coordinates introduced in this paper and to see, in particular, which linear combinations 
become irrelevant. 

We have also not pursued the study of other models, the Potts model, the n-vector model 
and so on. The examples of Section |^ indicate a surprising sensitivity to the definitions that 
it would be useful to examine further. For the Ising model we made, more by good luck than 
good management, a particularly happy choice which it is not utterly clear how to generalize 
to other models. 





(57) 



Kq{tp2,'ipi)dfi{ip2)dfi{ipi) = dmq{i)2,ipi) 
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Appendix. 



The present work contains simulations of both quahtative and quantitative nature. We 
aimed in most of Sections to || to provide numerical results reproducible to the precision 
of statistical errors. It is therefore important that we be precise about our conventions. 
Though many details are given in the text we complete them here with technical additions. 
In Paragraphs p.4| and O and in Section |^ the work is mostly qualitative and the reader 



who wants to examine these matters further will need to devise his own experiments. 
Distribution ttid- 

According to the principles of the introduction, each possible function h lying above a given a 
is to be assigned the same measure. This principle has to be incorporated into the programs 
locally. For example, there are two possibilities for the configuration of jump lines (or level 
curves) passing through the center of the configuration appearing in the first row of Figure 
|56| . They are chosen with equal probability. Since the curves are constructed one at a time 
by adjoining edges, when we first adjoin an edge passing through the center we then turn 
to the left or right with equal probability 1/2. The next time we pass through the center 
there is no choice; there is only one unused successor remaining. For a triangular lattice, 
there are no ambiguous configurations. For the hexagonal lattice, all combinations of + and 
— around vertices of the dual lattice lead to at most two possible choices of jump lines (and 
they are then treated as in the square lattice) except for the configurations in the two last 
rows of Figure |56| for which there are five possible local configurations of jump lines. Each 
will then have the probability |. As a consequence, when a curve first passes through the 
center of this configuration it continues on a straight line with probability | (which then 
leads necessarily to one of the configurations in the bottom row) or makes a sharp reverse 
turn to the left or the right with equal probabilities |. If the first curve through the center is 
straight, the following curves are determined. Otherwise the next curve, which may very well 
be a continuation of the first, returning after perhaps extensive wandering, has two options, 
each chosen with probability |. 



For the square lattice, two examples of the random determination of h occur in Figure |57 
where a configuration was drawn together with the jump lines of h. If the site at the bottom- 
left corner has coordinate (1, 1), then four clusters meet at (9i,2|) and at (13|,5|). In the 
first occurence, the two minus-clusters are joined and, in the second, they are separated. By 
definition the jump lines occur on edges dual to lattice bonds. Their vertices were rounded 
in this figure to show clearly the difference between joining and separating. The jump lines 
that wrap around the cylinder are indicated by dashed lines. 

The restrictions of h (on both the cylinders and the squares) were taken along several 
curves C. For the square lattices the curves were taken along lines of sites so that the 
intersection with dual bonds is unambiguous. The triangular lattices were oriented such 
that longitudinal lines with sites had one site per mesh unit. The longitudinal lines with 
sites of the hexagonal lattices had the pattern site-site- vacant repeated over every three-mesh 
cycles. For the triangular and hexagonal lattices the conformal images of the curves Ci on 
the cylinder never contained dual bonds parallel to them. They were however moved slightly 
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Figure 57: A configuration on a 11 x 21 cylinder with the jump lines of h. 
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to the closest position where their intersections with dual bonds were equally spaced. For 
the curves C = Cq at the boundary they were chosen as the curve closest from the boundary 
satisfying the previous requirement. 

We also measured the distributions on a disk of radius r = 300.2. The center of the disk 
was a site. All the sites inside the disk of radius r, and only these, were thermalized. Some 
of the boundary sites had three neighbors, others only two. We then determined an effective 
radius r^f^ as the radius of the largest circle that intersects only dual bonds associated to 
sites in the disk. It turned out to be r^s = 299.50. The restriction of h at the boundary was 
obtained along the circle of radius rgg — e with e = 0.001. The jumps in H are of the form 
±7iS{9 — 9o) where 6*0 is the position of the intersection on the circle of the dual bond with 
the curve C. The exact positions of all the intersections with dual bonds were determined 
and used to compute Fourier coefficients. The radius of inner circles were determined as 
fractions of {r^s — e). For example Cq and Ci are at 8 mesh units from one another on the 
397 X 793 cylinder and the radii of the corresponding curves on the disk should be [ves — e) 
and 0.8811(reff - e) since e'^'^-^/^^^ ~ 0.8811. 

Initial thermalization was provided by a few thousand Swendsen-Wang sweeps starting 
from a random configuration for the smaller cylinders, by 5000 sweeps for 397 x 793 and by 
10000 or more for 793 x 1585. After the initial thermalizations, measurements were taken 
every third Swendsen-Wang sweep for all the cylinders, except for the 793 x 1585 for which 
we used a 5-sweep cycle. A quick time-series analysis indicated that these cycle lengths 
insured proper statistical independence of consecutive measurements. The pseudo-random 



number generator was the one proposed by Tezuka and L'Ecuyer in [p^L 



Programs for the square lattice on the cylinders and on the disks were written indepen- 
dently by at least two of the authors and errors were chased down until measurements agreed 
within the statistical errors reported in the text. 

It might not be clear, on reading the main text, why certain data are given for some of the 
geometries studied and not for others. It is because the results for the various runs were kept 
in two different formats. For the first format, the observed values of each random variable 
or were grouped into 401 bins of equal size. Thus for each random variable, 401 
nonnegative integers were stored. The width of the bins used during the first months was 
a little too narrow and some of the values fell outside the range covered. Later, in the final 
months, the bin width was adjusted to avoid this difficulty. When we used the first format, 
we also kept, most of the time, the sum, the sum of the squares and the sum of mixed 
products of the random variables. These allowed us to calculate accurately the two first 
moments of the distribution and the correlation coefficients. The second format was more 
thorough. For each configuration of spins a, we constructed one possible h and recorded 
its restriction to the various curves C, not simply to the extremities of the cylinder or its 
median. When we realized that the conformal invariance might hold in the sense of Section 
0, we kept the restriction of h to more curves. These data bases, with only the boundary as 
C or with several curves C for each configuration, are sizeable even when compressed (a few 
Gigabytes in all). We generated one for the cylinder 397 x 793 with eleven curves C. For 
59 X 401 and 157 x 1067 we only kept the restriction of h to one extremity. With the second 
format it was possible to test various assertions that we could not have anticipated without 
the experience gained from the experiments, but the first required far less memory, so that 
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more sizes were examined. 



Crossings. 

To determine the aspect ratio of the rectangles where the crossings are measured we have 
used the width and the height of the smallest rectangle that contains the sites of the lattice 
considered. For example, for the rectangles of LV x LH sites of the square lattice, the ratio 
is LH/LV. We used here the orientations of the lattices used for the measurement of m^i 
(see above). 

For Tih, Tiv and Hhv on rectangles, crossings started on + spins from one boundary and 
ended on + spins on the other. For and n:^, the crossings were required to reach the 
central meridian if it contained sites or, if it did not, to reach the line of sites just before. 
The dimensions of the rectangles for the square lattice were the same as those used for 
percolation crossings in ||LPPS|| . The results for the triangular lattices are given in Table 
VII. The dimensions {LV, LH} that were used are 

{586,69}, {596,73}, {566,73}, {582,79}, {540,77}, {554,83}, 

{566,89}, {616,103}, {508,89}, {474,87}, {504,97}, {526,107}, 

{456,97}, {476,107}, {438,103}, {440,109}, {424,111}, {410,113}, 

{512,147}, {420,127}, {400,127}, {400,133}, {354,125}, {348,129}, 

{398,155}, {336,137}, {314,135}, {386,175}, {318,151}, {302,151}, 

{310,163}, {276,153}, {366,213}, {270,165}, {310,199}, {356,239}, 

{258,183}, {248,185}, {282,221}, {324,267}, {232,201}, {232,211}, 

{210,201}, {200,201}, {224,237}, {196,219}, {196,229}, {190,233}, 

{184,237}, {184,249}, {176,251}, {288,431}, {164,259}, {176,291}, 

{152,265}, {156,287}, {148,285}, {210,425}, {152,323}, {278,625}, 

{152,359}, {148,367}, {132,345}, {140,381}, {126,361}, {182,551}, 

{132,421}, {110,367}, {116,409}, {112,413}, {126,491}, {106,433}, 

{110,471}, {154,691}, {108,515}, {112,561}, {96,505}, {116,641}, 

{92,535}, {108,661}, {90,573}. 

The ratio r is given by r = 2LH/^/3LV as LH and LV count the number of lines and 
columns of sites. 

For TT/i, 7r„ and iThv on the disk, crossings started from and ended on sites in the annulus 
between r = 300.2 and r — ^/2. The crossings for vr^ and vr^ had to reach the central 
diameter. 

On cylinders the crossings between the curves Ci started from and ended on the curves. 
On the disk the five curves were chosen at radii f, 0.8811f, 0.7763f, 0.6026f, 0.3632f with 
f = 300.2 — ^/2. The crossings from Ci to Cj (r^ > rj) started outside the outer curve Q 
and ended inside the inner Cj. 

The programs for all lattices and geometries were written by two of us and checked until 
they agreed within the statistical errors for a sample larger than 10^ even though most 
crossings were measured with samples of ~ 200K. (See Section ^ for the samples used for 
the various lattices and geometries.) 
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The phase x. 



The phase x measured by the experiments is described in Section ^. For Figure |2^, results 
from cyhnders of the following sizes were plotted 

{59,27}, {59,37}, {59,47}, {59,61}, {59,73}, {59,93}, 
{59,119}, {59,147}, {59,179}, {59,211}, {59,249}, {59,283}, 

all with at least 400K configurations each, and 

{117,53}, {117,73}, {117,95}, {117,123}, {117,145}, {117,187}, 
{117,239}, {117,293}, {117,357}, {117,421}, {117,499}, {117,565}, 

with at least 600K configurations. 

The distribution of the random variable x is also used to obtain the ratios b/a through 
the constrained integrals ( ^I]) and (|3^). The errors on the ratios b/a appearing in Table 



VIII are difficult to evaluate as the numbers a and b are the local maxima of a smoothed 
distribution. For the integral (0), the most difficult to measure, the samples varied between 
31K and 50K. After experimentation with various smoothing parameters we think that the 
two first digits of the ratios b/a for the case constrained/ constrained are exact. The accuracy 
for the other cases is far better, the samples being at least 85K for the constrained/fixed and 
300K for the fixed/fixed. 



The correlation {{h{p) — /i(0))^). 

These correlations can be measured in a straightforward way using the above definition of h 
and the details in the text. 
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